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axis-switching  of  the  rings  )s  involved  in  the  process  is  seen  from  the  complex  topology  of  the  jet  shown  in  figure  (a). 
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Computational  Fluid 
Dynamics  - 
Navy  Perspective 

Spiro  G.  Lekoudis,  Guest  Editor 
Office  of  Naval  Research 


Abstract 

Computational  Fluid  Dynamics  (CFD)  combines  numeri¬ 
cal  algorithms,  fluid  dynamics,  and  computer  software  and 
hardware,  with  the  goal  of  predicting  flows.  With  rapid  devel¬ 
opments  during  the  last  three  decades,  CFD  changed/is  chang¬ 
ing  the  design  process  for  aircraft,  ships,  submersibles  and 
weapons.  CFD  will  be  important  to  the  Navy,  even  if  it  does 
not  impact  future  developments  as  prominently  as  it  did  in  the 
past.  This  is  the  case  because,  currently,  there  are  no  other 
paths  to  efficient  design  of  Naval  systems  whose  performance 
and  stealth  are  significantly  influenced  by  flows. 

Introduction 

About  two  decades  ago  I  was  learning  Computational  Fluid 
Dynamics  (CFD)  on  the  job,  while  working  as  aresearch  engineer 
in  the  defense  industry.  No  university  curricula  included  CFD  as 
a  topic  in  the  early  1970s,  and,  as  is  the  case  with  novel  ideas,  the 
science  &  technology  community  was  not  sure  about  its  potential. 
I  was  importing  CFD  software  from  the  “best  available”  sources 
in  government  and  universities,  and  I  was  applying  them  to  my 
problems.  It  is  interesting  that  this  practice,  a  generation  later,  is 
still  not  universally  accepted  as  the  most  cost  efficient.  I  had  both 
time-constrains  as  well  as  technological  barriers  and  they  were 
as  important  as  the  need  to  publish  my  work.  I  had  to  make  sure 
that  design  engineers  were  pleased  with  the  capabilities  that  I  was 


offering  them!  Since  time  was  critical  and  such  capabilities 
were  in  their  infancy,  I  had  to  learn  to  access  the  “best  and  the 
brightest”  and  their  software  in  other  organizations.  CFD  is 
the  product  of  this  generation  and  it  is  interesting  to  contem¬ 
plate  what  CFD  will  offer  to  the  Navy  as  the  “CFD  generation” 
is  getting  closer  to  retirement!  ONR,  always  involved  with  the 
forefront  of  science  &  technology,  supported  members  of  this 
“CFD  generation”  early  in  their  careers. 

About  twenty  five  years  ago,  researchers  developed  the 
ideas  needed  to  compute  steady  subsonic  flows  with  embed¬ 
ded  supersonic  regions  using  the  assumption  of  potential  flow. 
In  less  than  a  decade  the  computation  of  such  flows  became 
common  in  the  industry,  and  CFD  changed  forever  the  way 
aircraft  are  designed.  Such  change  is  occurring  in  the,  much 
smaller,  shipbuilding  industry.  Although  empirical  data  domi¬ 
nate  the  shaping  of  marine  vehicles  several  exceptions  exist, 
such  as  the  use  of  CFD  in  designing  the  winged-keel  of  racing 
yachts.  The  change  in  the  hydrodynamic  design  process  is 
influenced  by  increased  technological  demands  that  are  not 
easily  met  by  the  “cut-and-try”  methods,  by  the  power  of 
computers,  and  by  the  emerging  ability  to  deal  with  vortical 
flows.  Unlike  the  cruise  condition  of  a  commercial  jet,  where 
very  significant  quantities,  such  as  lift,  can  be  approximated 
by  potential  flow  CFD,  current  supercomputers  as  well  as 
powerful  workstations  allow  efforts  at  simulating  complex 
flows,  such  as  interacting  streamwise  vortices,  crucial  to  Na¬ 
val  technology. 
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CFD  is  a  combination  of  fluid  mechanics,  numerical 
algorithms,  and  computer  software  and  hardware,  that  at¬ 
tempts  to  simulate  fluid  flows.  The  relative  importance  of 
these  three  components  of  CFD  has  been  changing  over  the 
years.  Algorithm  development  dominated  CFD  in  the  1970s, 
but  both  its  pace  and  impact  is  slowing  down,  as  the  parallel 
supercomputer  industry  enters  a  more  “mature”  stage,  and 
load  balancing  becomes  the  primary  concern  of  users  of  such 
machines.  The  emergence  and  availability  of  parallel  super¬ 
computers  in  the  1980s  made  them  a  major  driver  in  CFD 
development.  In  the  1990s  powerful  workstations  provide  new 
access  to  processing  power.  Algorithm  development  of  a  dif¬ 
ferent  kind  is  occurring.  These  new  algorithms  are  not  targeted 
at  flow  solvers  but  implement  ideas  from  optimal  control 
theory  to  design  using  CFD  and  require  solutions  of  the 
adgoint  Navier-Stokes  equations.  Such  ideas  deal  with  the 
inverse  problem  in  vortical  flows.  The  solution  to  the  inverse 
problem  produces  solid  object  shapes  with  prescribed  con¬ 
strains  on  flow  loading  distributions. 

A  roadblock  that  slows  the  ever  increasing  applicability 
of  CFD  is  our  inability  to  simulate  turbulence  at  high  Reynolds 
numbers  and  for  the  interesting  complex  geometries.  Thus,  the 
importance  of  understanding  the  fluid  dynamics  using  theory  and 
measurements  cannot  be  overemphasized,  and’all  three  compo¬ 
nents  are  crucial  for  further  developments.  Unlike  the  early  stages 
of  CFD  development,  where  the  emphasis  was  in  obtaining 
quantitative  information  for  engineering  purposes,  CFD  is  also 
being  used  to  obtain  scientific  information,  for  example,  guiding 
measurements. 

The  design  and  operation  of  surface  ships  and  underwater 
vehicles  require  understanding  of  phenomena  that  involve 
incompressible  flows,  two-phase  flows,  cavitation,  hydroa¬ 
coustics,  and  free-surfaces.  The  presence  of  moving  solid 
surfaces  adds  further  complication.  Predicting  and  controlling 
such  phenomena  can  be  crucial  to  the  efficiency  of  various 
systems  and  to  their  stealth.  In  general  the  Reynolds  numbers 
are  very  large,  and  CFD  is  tasked  to  link  subscale  laboratory 
data  with  “full  scale”  flows  around  real  vehicles.  The  Reynolds 
number  is  proportional  to  the  characteristic  velocity  and 
length,  while  the  Froude  number  is  proportional  to  the  charac¬ 
teristic  velocity,  divided  by  the  square  root  of  the  characteristic 
length.  Since  the  Reynolds  number  and  the  Froude  number 
scale  differently,  either  empirical  formulas  or  CFD  has  to  be 
used  for  “scaling  up”.  In  this  way  low  Reynolds  number  wave 
tank  data  are  used  to  infer  data  at  full  scale.  It  is  difficult  to 
achieve  confidence  in  this  process,  but  this  extrapolation  is 
probably  the  most  useful  function  of  CFD  in  ship  design.  CFD 
needs  to  provide  information  that  reduces  design  uncertainties, 
and  reduces  the  amount  of  testing  and  trial  and  error  proce¬ 
dures  that  are  associated  with  design.  There  is  another  role  for 
CFD  in  hydrodynamic  design.  CFD  can  be  used  as  the  extrapo¬ 
lation  tool  that  allows  exploration  or  radical  departures  from 
existing  designs. 


The  Navy  involvement  with  CPD  started  at  the  Office  of 
Naval  Research  (ONR)  and  at  the  Naval  Research  Laboratory 
(NRL).  ONR  supported  both  fluid  dynamics  and  numerical 
algorithms.  Several  successful  CFD  methods  were  supported 
by  the  Navy  at  the  early  stages  of  their  development,  before 
they  became  practice  in  the  US  and  abroad.  Examples  include 
panel  methods,  spectral  methods  and  full-potential  transonics 
in  the  1960s  and  early  1970s,  and  Euler  solvers,  multi-grid  for 
Navier-Stokes,  spectral  elements,  unstructured  grids,  and  CFD 
for  parallel  computers  later  on.  Research  in  Computational 
Physics  at  the  Naval  Research  Laboratory  focused,  amongst 
other  issues,  on  shock  waves  and  reacting  flows.  The  Naval 
Air  Warfare  Center,  the  Naval  Surface  Warfare  Center,  and  the 
Naval  Undersea  Warfare  Center  develop  and  implement  CFD 
software  for  various  Naval  technology  uses,  for  example  plat¬ 
form/weapon  loads  during  weapons  release,  both  in  air  and 
underwater. 

The  remainder  of  the  article  emphasizes  CFD  topics  that 
are  important  and/or  unique  to  the  Navy’s  mission.  The  Navy 
does  benefit  from  CFD  developments  in  general,  but  there  are 
areas  that  are  of  particular  importance.  Moreover  the  DoD 
emphasis  on  service  reliance  and  the  significant  CFD  activities 
at  NASA  and  industry  help  to  increase  the  focus  of  Navy 
activities.  There  is  a  variety  of  physics  that  has  been  incorpo¬ 
rated  into  CFD  software.  Panel  methods  and  boundary  layer 
methods  continue  to  be  refined,  used,  and  contribute  to  tech¬ 
nology,  but  will  not  discussed  because  they  are  mature  tech¬ 
nology.  The  focus  will  be  on  CFD  for  the  Navier-Stokes 
equations,  which  are  used  at  some  point  in  vehicle  design 
process  by  the  sophisticated  design  teams. 

Some  Challenging  Issues  for 
Navy  Relevant  CFD. 

The  incompressibility  condition  introduces  an  elliptic 
constraint  in  the  unsteady  (parabolic)  Navier-Stokes  equa¬ 
tions.  This  constraint  has  been  numerically  satisfied  in  several 
ways.  In  formulations  that  involve  vorticity,  the  pressure  is 
eliminated  and  the  elliptic  effect  is  maintained  via  auxiliary 
Poisson  equations,  or  by  other  means.  In  methods  that  involve 
primitive  variables  a  Poisson  equation  for  the  pressure,  or  a 
pressure-like  variable,  is  solved.  Early  formulations  date  back 
to  the  1960s,  although  Navier-Stokes  solvers  were  then  a 
matter  of  academic  curiosity.  Another  approach  is  to  mimic  a 
fully  parabolic  problem  by  introducing  artificial  compressibil¬ 
ity.  If  we  confine  our  attention  to  numerical  algorithms  for 
single  processor  computers,  compressible  (and  single  phase 
and  species)  flow  CFD  is  a  rapidly  maturing  subject.  Artificial 
compressibility  makes  some  of  the  compressible  CFD  meth¬ 
ods  candidates  for  the  incompressible  problem.  Using  com¬ 
pressible  flow  solvers  for  incompressible  flows  is,  in  general, 
not  efficient,  due  to  the  time  step  limitations  as  the  Mach 
number  approaches  zero. 
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Vortical  motions  and  turbulence  underneath  or  at  the  free 
surface  interact  with  the  free  surface,  generate  waves,  and  may 
include  complex  gas-liquid  flows.  The  chemical  effects  of 
surfactants  on  the  surface  tension  complicate  the  physics.  CFD 
developments  are  needed  to  model  such  flows,  both  underwa¬ 
ter  and  with  firee-surfaces  present.  Most  of  the  current  approxi¬ 
mations  used  involve  integral  parameters  for  the  interaction  of 
bubbles  with  the  surrounding  liquid,  and  the  modification  of 
wave  motions  from  turbulent  wakes.  Currently  parabolized 
Navier-Stokes  are  used  to  compute  the  evolution  of  ship 
wakes,  and  this  evolution  is  correlated  to  wave  generation.  The 
free  surface  motion  couples  with  the  dynamics  of  ship  mo¬ 
tions.  CFD  could  become  crucial  in  understanding  and  engi¬ 
neering  for  such  motions,  especially  at  high  sea  states.  Such 
problems  are  currently  beyond  the  capability  of  Navier-Stokes 
to  predict.  To  do  that,  the  sophistication  of  CFD  predictions  of 
large  free-surface  deformations,  including  turbulence,  need 
significant  improvements. 

Surface  ships  generate  rather  complex  flowfields.  Design 
considerations  involve  powering  requirements,  sea-keeping 
characteristics,  and  acoustic  and  nonacoustic  signatures.  Tech¬ 
nologies  related  to  powering  and  sea-keeping  characteristics 
require  knowledge  of  the  near  ship  flows,  especially  their 
dominant  (bulk)  features.  These  flows  have  very  high 
Reynolds  numbers,  sometimes  in  the  billions.  Such  Reynolds 
numbers  bring  to  focus  the  issues  of  roughness  effects  and  wall 
functions  in  CFD.  Dominant  feature  of  near  ship  flows  are  the 
free  surface  and  the  separated  vortical  flows  surrounding  the 
hull.  Current  flow  solvers  address  such  flows  about  realistic 
hull  shapes  with  fully  nonlinear  boundary  conditions  at  the 
free  surface  and  with  “conventional”  turbulence  models. 
These  Reynolds-Averaged-Navier-Stokes  (RANS)  solvers 
employ  mostly  algebraic  turbulence  models,  and  the  total 
wave  drag  predicted  agrees  with  that  infeired  from  the  meas¬ 
urements.  Computational  ship  hydrodynamics  is  an  emerging 
area  and  is  discussed  in  the  article  by  Ed  Rood  in  this  volume. 

It  is  a  real  challenge  to  predict  similar  flows  either  while 
the  ship  is  maneuvering,  or  at  high  sea  states,  unless  the 
influence  of  breaking  waves  and  free-surface  turbulence  is 
minimal,  and  this  is  a  research  subject  for  Navier-Stokes  CFD. 
Current  results  for  such  flows  are  generated  using  unsteady 
potential  flow  methods  with  fully  nonlinear  free-surface  con¬ 
ditions.  Because  of  the  very  high  Reynolds  number  unsteady 
flows  of  interest,  a  RANS  solution,  where  all  scales  are  mod¬ 
eled,  needs  to  be  investigated.  However,  due  to  large  scale 
unsteadiness,  it  is  not  clear  at  this  time  if  a  RANS  approach, 
or  a  Very-Large-Eddy-Simulation  (VLES)  approach,  where 
most  of  the  scales  are  modeled,  is  a  superior  “predictor”.  A 
major  additional  consideration  for  the  computation  of  viscous 
unsteady  flows  is  the  severe  increase  in  computing  resources 
required  to  resolve  some  of  the  unsteady  phenomena  with  time 
accuracy.  If  the  physics  issues  related  to  the  predictive  capa¬ 
bility  of  the  software  are  resolved,  the  significant  new  power 


of  parallel  supercomputers  will  allow  routine  design  proce¬ 
dures  to  be  developed  for  viscous,  unsteady,  three-dimen¬ 
sional  flows  around  complex  geometries.  However  the  fidelity 
of  such  predictions  depends  on  models  for  turbulence,  maybe 
more  so  than  in  “steady”  flows.  Currendy,  predictions  of 
unsteady  loads  via  RANS  solutions  are  expensive,  not  always 
reliable,  and  thus  a  research  subject. 

The  features  of  small  scale  free-surface  waves  are  impor¬ 
tant  in  remote  sensing  applications.  Thus  free-surface  turbu¬ 
lence  and  the  effect  of  surfactants  play  a  significant  role  in  the 
physics  of  the  small  scale  features.  Simulations  have  been  used 
to  examine  the  effects  of  large  scale  vortical  motions  on  the 
characteristics  of  the  free  surface.  Direct  numerical  simula¬ 
tions  of  turbulent  open  channel  flow  produced  turbulence 
statistics  on  a  free  surface,  with  an  unperturbed  mean  location. 
The  behavior  of  air-sea  interface  is  part  of  a  larger  area  where 
numerical  fluid  dynamics  plays  an  important  role.  Oceans  are 
full  of  dynamic  events,  and  the  Navy  is  very  much  interested 
in  capabilities  that  allow  prediction  of  the  dynamic  behavior 
of  the  oceans  at  all  scales.  Temperature  and  density  variations, 
bottom  geometry,  and  air-sea  interactions  influence  such  pre¬ 
dictions.  Given  the  difficulty  of  obtaining  detailed  measured 
data  in  the  ocean,  the  role  of  simulations  may  be  crucial  in 
advancing  our  capabilities. 

In  underwater  vehicles  the  designer  deals  with  elongated 
bodies,  where  separated  flow  is  present,  and  where  the  propul- 
sor  operates  while  immersed  in  viscous  rotational  flow.  These 
flows  have  to  be  predicted  and  controlled.  In  general  they  are 
more  complex  than,  for  example,  attached  three-dimensional 
boundary  layer  flows.  Major  considerations  include  acoustic 
radiation,  ability  to  listen,  and  ability  to  absorb  impinging 
acoustic  waves.  These  considerations  involve  flow  over 
acoustic  sensors,  fluid  loads  that  lead  to  structural  vibration 
and  radiation,  and  flows  overcoatings  that  absorb  sound.  CFD 
is  tasked  to  help  design  with  acoustics  considerations.  RANS 
computation  for  complex  submarine  flows  task  the  most  pow¬ 
erful  of  supercomputers,  but  can,  in  principle,  provide  time- 
averaged  and  deterministic  (i.e.  non-turbulent)  unsteady 
pressure  and  shear  loads. 

In  the  submarine  world  acoustics  is  extremely  important, 
and  flow  details  that  do  not  contribute  much  to  loads  can 
become  significant.  Moreover  the  unsteady  nature  of  the  flow 
has  to  be  considered.  This  unsteadiness  ranges  from  the  large 
streamwise  vortices  originating  at  comer  flows  and  moving 
while  a  submarine  is  turning,  to  the  complex  pressure  foot¬ 
prints  of  turbulent  structures  within  boundary  layers  over 
acoustic  sensors.  The  unsteady  interaction  of  wakes  from 
control  surfaces  with  propulsors  involve  complex  three-di¬ 
mensional  phenomena  that  CFD  is  tasked  to  predict.  The 
difficulties  faced  by  CFD  are  illustrated  by  the  following 
example.  Assume  that  a  clever  designer  could  “engineer  out” 
much  of  the  shape  that  leads  to  flow  complexity.  Then  with 
simplifying  assumptions,  for  example  potential  flow  CFD,  he 
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could  get  satisfactory  steady  and  (maybe)  unsteady  loads.  This 
would  be  only  a  good  first  step,  because  designing  for  acous¬ 
tics  would  require  additional  details.  These  details,  such  as 
surface  pressure  fluctuations,  necessitate  sophisticated  CFD 
models  that  can  deal  with  the  unsteady,  separated,  turbulent 
flows  around  modem  marine  propulsors.  RANS  seem  to  be 
the  lowest  order  CFD  model  that  could  give  accurate  distribu¬ 
tion  of  important  unsteady  flow  quantities  over  solid  surfaces. 
This  is  because  RANS,  with  an  appropriate  turbulence  model, 
can,  in  principle,  capture  flow  separation  from  smooth  sur¬ 
faces.  How  separation,  large  vortical  motions,  and  unsteady 
loads  contribute  to  noise  generation  by  solid  surfaces,  as  such 
surfaces  move  through  complex  flows.  An  appropriate  remark 
on  this  topic  is  that  we  should  not  confuse  our  inability  to 
model  turbulence  with  the  errors  associated  with  neglecting 
physics.  As  an  example,  poor  performance  by  a  turbulence 
model  can  be  addressed  in  a  systematic  way,  but  inadequate 
separation  criteria  from  ad-hoc  assumptions  require  new  as¬ 
sumptions.  Turbulence  and  modeling  are  discussed  in  the 
article  by  Pat  Purtell  in  this  volume. 

Hydroacoustic  design  capabilities,  currently  under  devel¬ 
opment,  would  benefit  from  CFD  improvements  that  can 
speed  up  the  computation  of  unsteady  flows.  Flow  loads  lead 
to  structural  vibrations  that  can  result  in  acoustic  radiation, 
which  invites  comparisons  with  aeroelasticity.  Aeroelasticity 
is  important  for  the  design  of  propellers  and  helicopter  blades 
and  involves  large  structural  deformations  that  may  require 
truly  interactive  procedures  for  prediction  and  control.  In  such 
procedures  the  structural  deformation  is  predicted  by  coupling 
the  equations  of  motion  of  the  deforming  structure  with  CFD. 
In  water,  the  structural  deformation  can  be  orders  of  magnitude 
smaller  to  that  in  air  by  comparison,  and  still  produce  un¬ 
wanted  acoustic  radiation.  Important  physics  for  acoustic  sen¬ 
sors  are  due  to  the  fact  that  turbulent  flow  events  at  high 
wavenumbers,  couple  with  similar  events  in  the  structure  that 
are  difficult  to  simulate.  Moreover  low  wavenumber  events  in 
the  structure,  that  are  easier  to  compute,  couple  with  low 
wavenumber,  and  much  lower  amplitude  events  in  the  flow, 
that  are  extremely  difficult  to  simulate  with  CFD.  Thus  this 
linear  flow/structure  interaction  problem  is  very  difficult  to 
analyze,  and  novel  computational  methods  could  contribute  to 
this  topic. 

Computational  hydroacoustics  is  in  an  early  stage  of 
development,  but  progressing  rapidly.  Future  coupling  of 
acoustic  “sources”  where  the  fluid  dynamics  and  the  structural 
vibrations  dominate,  with  propagating  and  scattering  fields, 
where  computational  acoustics  is  needed,  will  provide  the 
Navy  with  enhanced  understanding  and  novel  predictive 
methods.  Cavitation  and  two-phase  (bubbly)  flows  occur  in  a 
variety  of  situations.  The  acoustic  effects  of  cavitation  are 
important  to  the  stealth  of  Navy  systems,  hence  cavitation  is 
to  be  avoided.  It  relates  to  the  local  pressure  field,  which  in 
turn  depends  on  the  flow  near  solid  boundaries  and/or  wakes. 


The  development  of  bubbly  flows  can  be  influenced  by  the 
surrounding  turbulence  near  solid  boundaries.  Current  engi¬ 
neering  CFD  predictions  rely  on  empirical  mixing  of  informa¬ 
tion  known  from  single  phase  turbulence,  with  integral 
parameters  that  incorporate  the  density  variation  present  in 
bubbly  flows,  or  by  computing  “two  species”  flows.  This  is  an 
excellent  research  area  for  CFD. 

Both  surface  ships  and  submarines  have  to  be  able  to 
withstand  shock  waves  generated  by  a  hostile  weapon.  In  case 
of  underwater  explosions  both  the  shock  wave  and  the  sub¬ 
sequent  generation  of  a  large  gas  bubble  can  inflict  damage. 
Understanding  how  to  design  for  minimum  damage,  and/or  for 
maximizing  weapon  effectiveness  requires  CFD  coupled  with 
prediction  of  structural  deformation  and  damage.  This  and 
other  similar  nonlinear  flow/structure  interaction  problems  are 
very  difficult  to  analyze.  The  great  variety  of  the  time  and 
spatial  scales,  and  the  transient  physics,  increase  the  demands 
on  CFD.  For  example  simulating  an  explosion  near  a  free 
surface  involves  the  energy  releasing  bubble  gases,  the  bubble 
free  surface,  the  ocean  surface,  the  structural  modes  of  the 
ship,  and  the  structural  failure  criteria. 

Landing  on  an  aircraft  carrier  in  adverse  weather  condi¬ 
tions  requires  special  attention  and  knowledge  of  the  forces 
and  moments  applied  to  the  landing  vehicle.  CFD  is  tasked  to 
provide  information  about  such  forces  and  moments  by  com¬ 
puting  the  complex  flow  around  the  aircraft/carrier  deck  ge¬ 
ometry.  Computing  separated  vortical  flows  is  a  major  CFD 
thrust,  especially  as  we  begin  to  quantify  the  airflow  interac¬ 
tion  with  moving  control  surfaces  and  the  loads  generated 
from  weapons  release.  Modern  military  aircraft  were  influ¬ 
enced  by  CFD  technology,  and  they  are  examples  of  systems 
with  signature  requirements  as  design  drivers.  CFD  is  also 
impacting  the  evaluation  of  mechanical  and  thermal  loads  on 
turbomachinery.  Navy  Warfare  Centers  are  investigating  is¬ 
sues  related  to  missile  warfare  with  CFD,  such  as  vertical 
launch  systems  with  confined  plumes  and  airframe/missile 
interactions  during  weapons  release.  Missle  related  CFD  is 
discussed  by  Robert  Van  Dyken  in  this  publication. 

When  turbulence  dominates,  floW  control  is  not  only  a 
most  difficult  research  topic,  it  is  also  of  great  technological 
importance.  In  the  past  CFD  had  essentially  no  role  to  play  in 
flow  control  of  the  small  turbulent  scales.  The  advent  of  direct 
and  large-eddy  simulations  provide  new  opportunities  to  ex¬ 
amine  flow  control.  The  computer,  in  principle,  can  be  used  to 
speed  up  the  search  and  confine  the  parameter  space  to, 
hopefully,  an  area  of  technological  significance.  Turbulent 
shear  flows  are  easier  to  control  when  they  are  away  from 
walls,  as  opposed  to  turbulent  boundary  layers.  An  example  of 
an  effort  to  use  CFD  for  control,  at  the  research  stage,  is  CFD 
of  simulations  of  combustion  instabilities.  Such  simulations, 
together  with  measurements,  are  used  to  develop  active  com¬ 
bustion  control  strategies.  The  topic  of  computational  combus¬ 
tion  is  discussed  by  Kailasanah  and  Roy  in  this  volume. 
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Such  efforts  point  to  how  crucial  is  the  coupling  of  CFD 
with  measurements.  The  influence  of  small  turbulent  scales  on 
the  large  vortical  motions  at  high  Reynolds  numbers  is  a 
complex  issue,  and  improved  understanding  via  both  simula¬ 
tions  and  measurements  may  lead  to  technologically  important 
control  strategies.  Direct  Numerical  Simulations  of  turbulence 
(DNS),  although  have  solved  the  “turbulence  problem”  at  low 
Reynolds  numbers,  will  nor  become  engineering  practice 
(high  Reynolds  numbers)  in  the  foreseeable  future.  However 
DNS  may  help  with  novel  turbulence  control  approaches, 
where  local  “event  controlling”  schemes  may  be  explored.  The 
rapid  development  of  microelectromechanical  systems 
(MEMS),  and  neural  net  controllers  may  help  the  implemen¬ 
tation  of  such  approaches.  However,  several  turbulence  con¬ 
trol  approaches  looked  good  in  the  Lab  but  never  became 
practice  because  the  scaling  laws  were  not  in  their  favor . 

Additional  Remarks. 

All  of  the  topical  areas  that  were  discussed  depend  on  all 
components  of  CFD,  fluid  dynamics,  algorithms  and  soft¬ 
ware/hardware.  The  emergence  of  parallel  supercomputers 
creates  opportunities  to  develop  software  with  increased  per¬ 
formance.  Although  the  time  scales  for  developing  and  imple¬ 
menting  such  software  for  engineering  use  depend  on  the 
objectives  of  the  individual  organizations  involved,  there  is  a 
clear  path  to  vastly  increased  computational  power.  Another 
change  in  computer  hardware  will  also  affect  the  level  of 
sophistication  that  is  acceptable  for  design.  New  powerful 
workstations  are  increasing  the  cost  effectiveness  of  sophisti¬ 
cated  software. 

However,  as  we  demand  increased  capabilities  for  pre¬ 
dicting  complex  flows,  our  dependence  on  modeling  turbu¬ 
lence  is  increasing.  Most  of  the  remaining  challenges  in  flow 
prediction  cannot  be  overcome  unless  we  understand  turbu¬ 
lence  better.  The  fact  is  that  low  Reynolds  number  turbulence 
can  be  directly  simulated  today,  and  that  numerical  fluid 
dynamics  is  providing  information  not  possible  to  obtain  by 
measurements.  The  practitioners  and  the  sponsors  of  CFD 
work  should  make  special  efforts  to  identify  the  sources  of 
disagreement  between  computations  and  measurements.  For 
example,  there  no  longer  excuse  for  not  knowing  that  the 
computation  of  two-dimensional  flow  suffers  from  inadequate 
numerics,  including  lack  of  grid  resolution.  Then  the  errors 
associated  with  modeling  turbulence  are  easier  to  assess.  Tur¬ 
bulence  remains  probably  as  the  most  difficult  unsolved  prob¬ 
lem  in  the  mechanics  of  continuous  media.  We  have  better 
numerics  and  larger  computers  than  the  days  when  predicting 
attached  transonic  flows  over  high  aspect  ratio  transport  wings 
was  the  forefront  of  CFD.  However  in  those  days  the  accuracy 
of  the  turbulence  model  was  less  important;  we  were  primarily 
interested  in  the  evaluation  of  integral  quantities,  such  as  wing 
lift,  while  now  we  are  seeking  prediction  of  local  phenomena. 


such  as  heat  transfer  on  turbine  blades  and  species  concentra¬ 
tion  in  missile  plumes. 

The  development  of  appropriate,  predictive  CFD  meth¬ 
ods  for  high  Reynolds  number  turbulent  flows  is  the  central 
issue  of  CFD.  Coupled  with  the  interdisciplinary  missions  that 
modem  CFD  is  tasked  with  brings  two  major  barriers  into 
focus.  These  barriers  cannot  be  overcome  by  better  numerical 
algorithms  and  more  powerful  parallel  supercomputers  alone. 
One  barrier  is  due  to  lack  of  good  models  for  turbulence,  and 
for  associated  multidisciplinary  physics,  for  example,  turbu¬ 
lent  combustion.  Theory  is  crucial  in  accomplishing  this,  and 
the  importance  of  new  theoretical  developments  cannot  be 
overemphasized.  The  second  barrier  is  due  to  lack  of  measure¬ 
ments  at  high  Reynolds  numbers,  so  that  the  correct  physics 
can  assist  the  modeling  efforts.  Facilities  such  as  the  ARPA 
“Superpipe”  at  Princeton  and  the  new  NASA/Industry  high 
Reynolds  number  facilities  (at  the  drawing  board)  will  con¬ 
tribute  to  obtaining  such  data.  It  is  CFD  validation,  the  process 
by  which  software  and  measurements  are  critically  examined 
to  establish  the  capabilities  and  limits  of  the  software  that  will 
promote  further  acceptance  of  CFD.  Such  validation  efforts, 
costly  in  resources  and  time,  are  absolutely  necessary. 

Even  if  we  only  partially  succeed  in  uanslating  new 
knowledge  about  turbulence  to  engineering  tools,  CFD  has 
already  contributed  in  another  way.  CFD  driven  knowledge 
about  numerical  algorithms  and  supercomputer  use  is  being 
incorporated  into  other  scientific  and  technological  areas,  for 
example,  computational  electromagnetics.  CFD  matured  al¬ 
gorithms  that  are  used  to  improve  manufacturing  processes, 
such  as  the  design  and  analyses  of  computer  chips,  to  under¬ 
stand  the  behavior  of  structures  under  impact  loading,  and  to 
improve  combustion  and  energy  conversion  related  technolo¬ 
gies.  These  few  examples  point  to  current  and  future  activities 
that  will  increase  the  Navy’s  ability  to  perform  its  mission. 
Although  such  activities  are  not  within  my  definition  of  CFD 
they  benefit  from  it,  because  CFD  stands  out  as  an  example  of 
how  large  computers  proved  useful  to  DoD  in  important 
product  development. 


Summaiy. 

Computational  Fluid  Dynamics  changed/is  changing  the 
platform  design  process  and  is  influencing  several  other  Navy 
relevant  technologies.  This  occurred  not  only  because  of  the 
capabilities  of  new  supercomputers,  but  also  because  CFD 
reduced  the  cost  of  obtaining  technologically  significant  quan¬ 
tities.  However,  as  is  the  case  with  all  areas  of  science  and 
technology  that  compete  for  resources,  its  future  significance 
depends  on  continuing  technological  impact.  There  are  major 
obstacles  in  CFDs  very  successful  path.  Current,  planned  and 
envisioned  computer  hardware  improvements  can  help  to 
overcome  these  obstacles,  but,  by  themselves,  cannot  remove 
them.  Such  obstacles  are  primarily  due  to  lack  of  appropriate 
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models  and  of  appropriate  measurements  for  high  Reynolds 
number  turbulence,  and  for  the  problems  generated  when  such 
turbulence  affects  other  phenomena. 

Both  the  engineering  and  the  scientific  questions  that 
challenge  CFD  today  are  very  demanding.  In  the  future,  CFD 
will  be  part  of  simulation  science  and  simulation  science  will 
be  dominant  in  concurrent  engineering  and  simulation  based 
design  (SBD).  The  vision  of  combining  fluid  dynamics  with 
numerical  algorithms  and  powerful  supercomputers  has  al¬ 
ready  helped  the  design  of  new  platforms.  Moreover,  there  are 
currently  no  alternatives  to  CFD  that  satisfy  Navy  needs  in  the 
topics  discussed  in  this  article. 
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Abstract 

The  field  of  computational  ship  hydrodynamics  has  re¬ 
cently  emerged  as  a  result  of  the  application  of  contemporary 
computational  fluid  dynamics  to  ships  and  submarines.  Com¬ 
putational  ship  hydrodynamics  incorporates  a  validated  nu¬ 
merical  process  that  models  as  one  system  the  interacting 
hydrodynamic  phenomena  in  the  flow  around  a  marine  vehi¬ 
cle.  Hydrodynamic  flows  around  marine  vehicles  are  unique 
because  the  vehicles  are  characteristically  drag-producing  and 
become  dynamic-lift-producing  only  in  turns  or  at  exception¬ 
ally  high  speeds.  The  vehicles  have  complex  responses  to 
control  surface  loads  as  a  consequence.  Often  ship  hydrody¬ 
namic  flows  are  bounded  by  free-surfaces  and  require  a  solu¬ 
tion  for  the  boundary  deformation  as  well  as  the  field  variables 
such  as  the  velocity  and  pressure.  Usually  marine  propulsors 
are  located  in  the  wake  of  the  hull  and  appendages,  as  well  as 
the  wake  of  their  own  shafts  and  struts.  There  is  a  strong 
interaction  between  the  flow  around  the  hull  and  the  flow 
through  the  propulsor  and  the  rudder.  As  in  many  other  engi¬ 
neering  areas,  the  availability  and  the  exploitation  of  super¬ 
computer  resources  has  provided  the  needed  breakthrough  to 
rationally  analyze  and  design  complex  ship  hydrodynamic 
systems.  Recent  computational  milestones  include  the  steady 
hydrodynamic  flow  around  propelled  surface  ships  and  the 
design  of  complex  integrated  stem/propulsor  systems.  Nu¬ 
merical  capability  has  produced  results  that  are  more  detailed 
than  existing  measurement  data;  hence  there  is  a  need  for  new, 
more  accurate,  highly  resolved  measurements  of  the  unsteady 
flow  around  large  model  ships  in  towing  tanks  and  maneuver¬ 


ing  basins  in  order  to  validate  and  calibrate  the  computations. 
In  the  near  future  computational  ship  hydrodynamics  will  be 
an  integral  part  along  with  traditional  towing  tank  tests  for 
design  concepdon  and  evaluation.  In  the  not-too-distant  future 
the  naval  architect  can  expect  to  use  a  numerical  towing  tank 
as  the  primary  design  tool,  with  the  towing  tank  providing 
calibradon  and  limited  validadon  at  model  scale  for  the  final 
design. 

Introduction 

Fluid  dynamics  is  a  venerable  subject  occupying  the 
scientific  minds  of  mathematicians,  physicists,  and  engineers. 
Fluid  dynamics  as  a  modem  science  arguably  dates  from  1500 
with  Leonardo  da  Vinci’s  flow  visualizations  (Lugt,  1983). 
Four  hundred  years  later,  the  modem  age  of  the  supercomputer 
has  been  a  revolutionary  period  for  fluid  dynamics.  As  recent 
as  1980,  only  visionaries  anticipated  that  numerical  simula¬ 
tions  of  the  physics  of  such  flows  as  three-dimensional  turbu¬ 
lent  jets  interacting  with  the  free  surface  (Mangiavacchi,  et  al., 
1994)  and  wave  interactions  with  turbulent  patches  on  the  free 
surface  (Dommermuth,  et  al.,  1994)  would  be  commonplace 
in  1995.  This  article  is  about  how  computational  fluid  dynam¬ 
ics  is  impacting  the  field  of  naval  ship  hydrodynamics,  and 
about  how  and  why  the  emergence  of  computational  ship 
hydrodynamics  from  a  numerical  supplement  to  a  numerical 
complement  to  towing  tank  tests  will  lead  to  a  tum-of-the-cen- 
tury  revolution  in  the  design  process. 

Ship  hydrodynamics  is  similar  to  airplane  aerodynamics 
in  many  respects,  especially  in  the  needs  to  properly  grid  the 
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Figure  1. 

Characteristic  Free-Surface  Features  of  the  Hydrodynamic 
Flow  Around  a  Ship. 


vehicle  geometry,  to  utilize  an  appropriate  turbulence  model, 
and  to  account  for  flow  separations,  and  unlike  airplane  fluid 
dynamics  in  many  other  respects.  For  instance,  ships  and 
submarines  are  essentially  drag  bodies  with  the  propulsor  in 
the  wake  of  the  vehicle  whereas  airplanes  are  lifting  bodies 
with  the  propulsors  in  front  of  the  body.  Flow  features  such  as 
the  free  surface  and  two-phase  flows,  shown  in  Figure  1 ,  must 
be  included  in  computational  ship  hydrodynamics  whereas 
those  features  are  absent  from  airplane  aerodynamics.  Marine 
propulsors  are  characterized  by  complex  geometry  featuring 
strongly  three-dimensional  flows,  whereas  aircraft  propulsors 
feature  flows  more  two-dimensional  in  character.  These  dis¬ 
similarities  between  ship  hydrodynamics  and  airplane  aerody¬ 
namics  have  led  to  the  specialized  field  of  computational  ship 
hydrodynamics. 

This  article  is  directed  at  scientists  and  technologists  with 
little  background  in  fluid  dynamics  as  well  as  to  computational 
ship  hydrodynamicists  interested  in  the  larger  role  of  numeri¬ 
cal  methods  in  science  and  technology  and  Navy  R&D.  The 
article  continues  the  perspective  generated  by  Lekoudis  (1993, 
1995)  but  focused  here  on  ship  hydrodynamics.  The  article  is 
long  on  the  role  and  technological  implications  of  computa¬ 
tional  ship  hydrodynamics  and  short  on  the  description  of 
numerics.  It  is  an  overview  with  selected  illustrative  examples 
of  recent  milestone  accomplishments  including  computations 
of  the  fully  nonlinear  free-surface  turbulent  flow  around  a 
combatant  surface  ship,  the  computation  of  the  propeller-hull 
interaction  and  in  particular  the  elusive  “effective  wake”,  and 
the  computation  of  the  flow  around  a  sailing  yacht  at  yaw 
including  the  lifting  surface  keel.  In  addition,  examples  of  very 
high  resolution  simulations  of  selected  phenomena  in  the  flow 
around  a  ship  are  presented.  These  high  resolution  simulations 


are  used  to  obtain  coefficients  for  models  embedded  in  the 
numerical  method  for  predicting  the  overall  flow  around  the 
ship. 

The  Equations 

The  future  for  computational  ship  hydrodynamics  is  in  the 
numerical  solution  of  the  unsteady  Reynolds-Averaged 
Navier-Stokes  (RANS)  equations  in  which  the  complete 
Navier-Stokes  equations  for  a  fluid  are  approximated  by  as¬ 
suming  the  flow  is  incompressible  and  by  modeling  the 
Reynolds  stresses  (the  “turbulence").  It  turns  out  that  these 
assumptions  do  not  practically  affect  the  validity  of  the  equa¬ 
tions  for  ship  hydrodynamics.  The  equations  are  solved  subject 
to  the  boundary  conditions  presented  by  the  ship  hull  and  the 
environment,  such  as  the  deformable  free  surface  and  the 
muddy  ocean  bottom. 

The  unsteady  RANS  equations  are  the  momentum  equa¬ 
tions,  or  the  equations  of  motion,  for  the  fluid: 

=  ^  +  V  <u'u'> 

dt  p 

where 

<u'u'>  =  UiUj  ( /  =  1 ,2,3;  j  =  1 .2,3  ) 

is  the  nine-component  velocity  correlation  tensor  for  the 
turbulent  velocity  fluctuations;  it  is  the  Reynolds  stress  ten¬ 
sor.  These  equations  (one  for  each  coordinate  direction) 
admit  important  flow  characteristics  including  stratifica¬ 
tion,  temperature  distribution,  and  bubble  clouds  with  the 
appropriate  coupled  equations. 

The  unsteady  RANS  equations  are  obtained  by  averaging 
over  characteristic  times  longer  than  the  turbulence  charac¬ 
teristic  time  yet  shorter  than  the  characteristic  time  for  ship 
motions.  The  equations  predict  motions  corresponding  to  ship 
turning,  propeller  rotation,  and  wave  encounter  while  account¬ 
ing  for  the  averaged  effect  of  small  scale  random  motions.  The 
Reynolds  stress  represents  the  average  effect  of  the  small 
scales,  and  is  obtained  by  averaging  the  nonlinear  velocity 
term  in  the  Navier-Stokes  equation.  Hence  it  is  an  inertial  term 
rather  than  a  viscous  term.  It  is  the  Reynolds  stress  for  which 
models  are  sought  in  the  quest  to  “solve  turbulence”  (Purtell, 
1995). 

The  RANS  equations  are  solved  subject  to  mass  conser¬ 
vation  and  the  boundary  conditions.  The  boundary  conditions 
are  specified  on  all  the  boundaries  of  the  fluid,  including  the 
ship  hull,  the  deformable  free  surface,  and  the  bottom  and  sides 
of  the  fluid  domain.  Often  the  sides  are  considered  to  be 
sufficiently  far  from  the  ship  that  they  do  not  have  a  direct 
effect  on  the  flow  around  the  ship.  In  this  case  the  boundary 
of  the  fluid  is  the  boundary  of  the  numerical  grid,  and  condi- 
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Figure  2. 

Photograph  of  the  Cotpplex  Stern,  Propeller,  Strut,  and  Rudder 
Geometry  Leading  to  Complicated  Three-Dimensional  Row  In¬ 
teractions. 
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tions  such  as  no  wave  reflection  and  no  streamwise  gradients 
of  flow  quantities  are  imposed  at  the  boundary  of  the  grid. 

The  hull  boundary  conditions  are  expressed  as  a  conse¬ 
quence  of  the  “no-slip”  condition  which  specifies  that  the  fluid 
at  the  surface  of  the  hull  assumes  the  velocity  of  the  hull.  This 
leads  to  regions  near  the  hull,  and  in  shed  vorticity,  in  which 
the  flow  exhibits  very  large  momentum  gradients.  There  is 
then  a  wide  range  of  length  scales  to  be  resolved  by  computa¬ 
tions  of  the  flow.  The  range  of  scales  is  effectively  defined  by 
the  Reynolds  number  (the  ratio  of  the  inertial  to  the  viscous 
forces),  with  the  range  increasing  rapidly  with  increasing 
Reynolds  number.  The  Reynolds  number  is  small  in  laboratory 
experiments,  moderate  in  towing  tank  tests,  and  very  large  for 
full  scale  ship  hydrodynamics.  These  high  gradient  regions 
cannot  be  ignored  because  they  are  the  boundaries  between 
the  bulk  fluid  field  and  the  hull  and  hence  transmit  the  fluid 
force  to  the  hull,  and  because  they  are  the  sources  for  various 
flow  interactions  as  the  high-gradient  fluid  is  advected  down¬ 
stream. 

The  free  surface  is  a  material  surface  boundary  to  the 
liquid  flow  domain  on  which  additional  boundary  conditions 
are  imposed.  Those  conditions  range  from  the  simplest,  in 
which  the  surface  is  considered  to  be  a  constant  pressure 
boundary  with  zero  viscous  stress,  to  more  complex  in  which 
the  wind-induced  viscous  stress  is  included  in  the  boundary 


condition.  Even  in  the  simplest  case,  the  presence  of  the  free 
surface  complicates  the  computation  because  the  location  of 
the  free  surface  must  be  determined  as  part  of  the  solution  to 
the  problem. 

Submarines  are  included  as  “ships”  in  this  article.  Al¬ 
though  in  many  cases,  such  as  deep  submergence,  free  surface 
effects  are  absent,  the  submarine  hydrodynamics  are  charac¬ 
terized  by  their  own  peculiar  features  attributable  to  the  wake 
inflow  to  the  propulsor  and  to  unusual  motions  caused  by  high 
incidence  angles  to  the  flow  and  irregular  flow  separations 
during  maneuvers. 

The  Essence  of  Computational 
Ship  Hydrodynamics 

Computational  ship  hydrodynamics  includes  the  collec¬ 
tion  of  numerical  methods,  the  verification  of  the  methods 
through  study  of  the  convergence  characteristics  of  the  nu¬ 
merical  methods,  and  the  validation  of  the  computed  results 
by  physical  measurements.  Convergence  studies  and  valida¬ 
tion  measurements  are  necessary  because  the  numerical  meth¬ 
ods  are  alternative  representations  of  the  equations  for  the 
flow,  and  feature  dynamics  dictated  by  their  method  of  solu¬ 
tion  as  well  as  by  the  flow  physics.  It  is  possible  for  a  method 
to  converge  to  the  wrong  answer,  or  to  give  a  seemingly 
reasonable  answer  which  is  grid  dependent  and  hence  not 
converged. 

The  naval  requirement  is  to  understand  the  details  of  the 
flow  around  marine  vehicles  at  several  levels.  At  the  coarsest 
level,  knowledge  of  the  accurate  time-dependent  evolution  of 
the  macroscale  flow  over  the  body,  especially  including  inter¬ 
ference  effects  of  appendages  and  propulsors,  is  required  to 
determine  forces  and  moments  during  maneuvers.  Even  at  this 
coarse  level  the  complexity  of  the  ship  geometry,  an  example 
of  which  is  shown  in  Figure  2,  leads  to  very  complicated 
three-dimensional  flow  interactions.  The  computed  forces  and 
moments  must  be  sufficiendy  accurate  that  they  can  be  used 
to  predict  and  control  surface  ship  motions  in  response  to  a 
seaway,  and  submarine  trajectories  for  extreme  maneuvers, 
such  as  might  be  expected  under  emergency  situations.  At  a 
more  detailed  level,  a  prediction  of  the  accurate  time-depend¬ 
ent  evolution  of  the  flow  details  into  the  propulsor  is  required 
for  acoustic  design.  These  details  include,  for  example,  the 
spatial  distribution  of  the  inflow  to  the  propeller.  This  spatial 
distribution,  which  appears  as  a  time-dependent  flow  in  the 
rotating  coordinate  frame  of  the  propeller,  is  required  to  pre¬ 
dict  vibrations  of  the  propeller.  Inherent  in  the  computation  of 
the  inflow  for  marine  vehicles  is  the  need  to  correctly  predict 
the  evolution  of  the  vortical  wakes  produced  by  upstream 
disturbances  such  as  appendages  and  even  by  thq  hull  if  it  is 
inclined  to  the  flow.  Hence  the  computational  method  must 
not  be  overly  dissipative.  At  the  most  detailed  level,  the 
accurate  simulation  of  the  microscale  evolution  of  the  flow  is 
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required.  Such  simulations  are  used  for  example  to  reproduce 
the  flow  as  a  numerical  measurement  for  the  calibration  of 
turbulence  models. 

The  process  of  computing  the  hydrodynamic  flow 
around  a  ship  begins  with  the  specification  of  the  ship 
geometry.  This  first  step  is  crucial.  The  ship  hull,  append¬ 
age,  and  propeller  geometry  defines  the  boundary  to  the 
flow;  small  changes  in  the  geometry  can,  under  appropriate 
circumstances,  lead  to  significant  changes  in  vehicle  per¬ 
formance.  Discontinuities  in  the  geometry  can  produce  dis¬ 
proportionate  disruptions  to  the  computed  flow.  The 
geometry  is  developed  initially  in  computer  aided  design 
applications,  and  includes  highly  resolved  interfaces  be¬ 
tween  such  components  as  the  hull  and  appendages. 

The  next  step  in  the  process  is  the  conversion  of  the 
flow  volume  and  boundary  surface  geometry  to  a  discrete 
grid.  In  this  step  the  continuous  ship  geometry  is  trans¬ 
formed  to  a  numerically  similar  array  of  points  in  space 
according  to  the  complexity  of  the  geometry  and  the 
anticipated  flow.  The  fluid  volume  and  the  deformable 
free  surface  are  also  transformed  to  a  three-dimensional 
array  of  points  in  space  representative  of  the  actual  continu¬ 
ous  fluid  domain.  The  creation  of  the  grid  is  a  major  con¬ 
sumer  of  resources,  especially  in  the  case  of  time-dependent 
flows  in  which  the  grid  must  be  adjusted  at  each  time  step 
to  account  for  the  evolving  flow,  such  as  the  changing  free 
surface  deformation  or  the  downstream  movement  of  vor¬ 
tices.  The  grid  must  be  sufficiently  fine  to  resolve  space  and 
time  details  but  must  facilitate  efficient  computations  and 
not  lead  to  unstable  or  nonconverged  solutions. 

In  the  next  step,  the  equations  of  motion  and  mass 
conservation  are  discretized  and  coupled  to  a  numerical 
solution  method.  The  discretization  and  solution  method 
can  be  any  one  of  a  number,  including  finite  volume  and 
finite  difference  methods.  The  discretization  is  a  numeri¬ 
cal  approximation  to  the  analytic  equations;  inversion  of 
the  discretized  equations  produces  a  mathematical  ex¬ 
pression  only  approximating  the  initial  equation.  It  is 
expected  that  the  solution  first  obtained  is  not  the  desired 
representation  of  the  actual  flow.  A  crucial  check  on  the 
approximation  is  convergence  with  changing  grid  size 
and  with  iterations  in  time.  The  numerical  solution  is 
“unstable”  if  the  computations  become  unrealistic  as  the 
convergence  is  pursued.  Stability  is  essential  for  the 
computations  to  produce  a  meaningful  answer  and  can 
sometimes  be  obtained  by  increasing  the  effective  damp¬ 
ing  in  the  numerical  solution  method.  A  stable  numerical 
procedure  is  “robust".  However,  there  is  a  critical  balance 
struck  between  the  damping  necessary  to  stabilize  the 
equation  solutions  and  the  marginal  stability  required  to 
retain  flow  details.  For  example,  robust  and  overly- 
damped  computations  will  cause  vortical  motions  to  de¬ 
cay  unnaturally  and  lead  to  false  conclusions  regarding  the 


effects  of  upstream-generated  vorticity  on  downstream  flow 
features. 

The  next  step  in  the  computational  process  is  the  presen¬ 
tation  and  communication  of  the  computed  results  as  useable 
information.  This  step  is  a  major  problem:  the  flows  are 
complicated  and  the  geometry  highly  complex,  and  the  data 
space  is  multidimensional  including  geometry,  velocity,  pres¬ 
sure,  temperature,  density,  and  derivatives  such  as  vorticity. 
Analysis  of  this  information  requires  very  sophisticated  visu¬ 
alization  processing,  that  can  render  iso-surface  contours  and 
perspective  views  “instantaneously.”  Such  “post  processing” 
is  always  accompanied  by  three-dimensioinal  visual  presenta¬ 
tion  of  the  data,  usually  generated  by  sophisticated  computer 
workstations.  Characteristically  such  information  appears  as 
visual  tracking  of  flow  particles  as  they  move  in  time  through 
color-coded  pressure  regions,  with  changes  in  views  and  per¬ 
spectives  being  made  “instantaneously".  Virtual  reality  will 
soon  have  an  impact  by  permitting  the  analyst/designer  to 
virtually  walk  around  the  flow  and  observe  phenomena.  This 
observation  is  crucial  for  the  full  exploitation  of  computational 
ship  hydrodynamics. 

The  final  step  in  the  process  of  computing  the  hydrody¬ 
namic  flow  around  a  ship  is  the  validation  of  the  results.  A 
limited  set  of  physical  measurements  is  required  to  provide 
assurance  that  the  compulations  are  useable  for  analysis  and 
design.  Within  a  reasonable  latitude,  such  measurements  can 
be  used  to  actually  calibrate  the  computed  solution;  for  exam¬ 
ple,  certain  turbulence  model  coefficients  may  be  set  by  the 
measurements.  The  process  that  led  to  the  validated  solution 
can  then  be  safely  repeated  to  develop  a  data  base  for  various 
conditions. 

In  time,  and  with  extensive  experience,  the  validated 
computational  ship  hydrodynamics  process  can  lead  to  a  cer¬ 
tified  process  able  to  be  relied  upon  for  absolute  engineering 
solutions.  Computational  ship  hydrodynamics  has  not  realized 
this  state  of  maturity. 

Towards  a  Numerical  Towing 
Tank 

Computational  fluid  dynamics  has  reached  a  level  of 
sophistication  and  validity  such  that  it  is  routinely  used  to 
compute  hydrodynamic  flows  around  various  components  of 
the  ship,  such  as  the  hull,  the  appendages,  and  the  propeller. 
The  next  important  step  is  to  combine  the  computational 
methods  into  one  or  more  comprehensive  computational 
methodologies  and  techniques,  or  computational  ship  hydro¬ 
dynamics.  Such  a  computational  ship  hydrodynamics  capabil¬ 
ity  is  technically  feasible  today  as  a  tool  useable  by  the  design 
community. 

Experience  with  successful  applications  of  computational 
fluid  dynamics  to  real-world  problems  has  demonstrated  the 
need  to  form  a  group  of  individual  specialists  dedicated  to 
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applications-driven  research  (Whitfield,  1995a).  This  man¬ 
agement  concept  is  contrasted  to  modified  matrix-manage¬ 
ment  often  advocated  but  usually  lacking  the  dedication 
characteristic  of  a  team.  In  particular,  experience  indicates 
that  the  Navy  is  best  served  by  a  hydrodynamics  center 
devoted  to  computational  ship  hydrodynamics.  This  center 
would  gain  experience  as  the  individual  specialists  (grid 
generation,  turbulence  models,  wake  simulation,  post-proc¬ 
essing,  etc.)  work  as  a  dedicated  team  on  a  variety  of 
problems.  The  Navy  has  established  the  Hydrodynam¬ 
ics/Hydroacoustics  Technology  Center  at  the  David  Taylor 
Model  Basin  in  Bethesda,  Maryland,  to  address  the  imple¬ 
mentation  of  computational  ship  hydrodynamics  in  the  de¬ 
sign  world. 

Of  all  the  hydrodynamic  components  on  a  ship,  the 
propeller  has  been  the  most  designed  by  numerical  compu¬ 
tational  methods.  Propeller  design  traditionally  includes  a 
substantial  reliance  on  computational  fluid  dynamics  in 
combination  with  experience.  Numerical  methods  ranging 
from  lifting  line  to  obtain  spanwise  blade  loading,  lifting 
surface  to  design  the  chordwise  loading,  and  (more  re¬ 
cently)  panel  methods  to  design  to  avoid  cavitation  are  used 
to  obtain  a  set  of  geometries  which  are  evaluated  in  water 
tunnels  for  applications.  Computational  ship  hydrodynam¬ 
ics  is  now  admitting  the  capability  to  predict  the  combined 
interaction  flow  for  a  propelled  ship.  It  is  anticipated  by 
naval  hydrodynamicists  that  the  role  of  computational  ship 


hydrodynamics  will  become  increasingly  more  visible  and 
credible.  The  rate  at  which  there  will  be  increased  utilization 
of  computational  ship  hydrodynamics  is  regulated  more  by 
perception  of  validity  than  by  technological  capability. 

Computational  ship  hydrodynamics  has  reached  a  level 
of  maturity  at  which  it  can  be  used  for  sorting  various  real- 
world  designs  at  least  for  a  relative  ranking.  The  expectation 
in  the  hydrodynamics  community  is  that  computational  ship 
hydrodynamics  can  and  should  be  a  major  tool  for  the  concep¬ 
tual  and  preliminary  designs  for  the  next  surface  combatant 
(SC-21),  and  that  this  tool  will  prove  essential  for  problem¬ 
solving  in  response  to  towing  tank  test  results  for  designs. 

The  computational  ship  hydrodynamics  approach  is  de¬ 
veloped  for  a  particular  range  of  geometry  and  operating 
conditions.  This  specialization  is  required  for  the  foreseeable 
future  because  of  the  sensitivity  of  the  results  from  computa¬ 
tional  methods  and  models  to  geometry  and  flow  conditions. 
The  methods  are  necessarily  marginally  stable  to  ensure  maxi¬ 
mum  fidelity  of  the  computations.  The  procedure  is  to  cali¬ 
brate  the  computations  with  a  limited  set  of  towing  tank 
measurements,  such  as  the  straight-ahead  operation  in  calm 
and  moderate  sinusoidal  waves  for  one  set  of  operating  con¬ 
ditions.  The  calibration  is  not  complete  because  the  model 
scale  flow  is  not  the  same  as  the  full  scale  flow,  in  part  because 
of  the  incongruence  of  modeling  both  the  gravitational  and 
viscous  forces  simultaneously.  Nevertheless,  the  model  scale 
flow  exhibits  features  similar  to  those  found  on  full  scale,  as 
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shown  in  Figure  3  comparing  full  scale  and  model  scale 
features. 

Once  validated  and  calibrated  for  the  design  condition, 
computational  ship  hydrodynamics  can  be  used  to  analyze  ge¬ 
ometry  modifications  and  off-design  operation  including  maneu¬ 
vers.  For  example,  wake-adapted  propulsors  designed  for 
minimum  unsteady  loading  in  the  straight  ahead  condition  may 
be  better  optimized  over  a  range  of  angles  of  incidence  produced 
by  maneuvCTS  or  by  waves.  Designs  can  incoiporate  the  range  of 
parameters  expected  under  normal  operation  rather  than  being 
focussed  on  the  “straight  ahead”  condition  with  heuristic  knowl¬ 
edge  incorporated  into  the  design. 

In  the  world  of  computational  ship  hydrodynamics,  there 
is  no  practical  limitation  on  the  scale  for  the  ship.  Therefore 
the  computed  model  scale  data  can  is  scaled  to  full  scale  values 
useful  for  design.  This  procedure  does  not  entirely  remove  the 
empirical  “correlation”  connecting  model  scale  and  full  scale 
results,  but  it  does  significantly  increase  the  confidence  in  the 
results  because  the  scaling  is  physics-based.  This  is  especially 
important  if  the  information  desired  lies  outside  the  historical 
parameter  range. 

Computational  ship  hydrodynamics  can  be  utilized  in 
several  ways  to  impact  design  and  analysis.  The  following 
several  examples  serve  to  illustrate  some  of  those  ways. 

Ship  motions,  and  in  particular  rolling  motions,  are 
greatly  affected  by  the  viscous  and  turbulent  damping  of  the 
motion  by  the  water.  This  damping  depends  on  the  Reynolds 
number,  which  is  very  large  for  full-scale  ships,  and  small  for 
model  scale  ships.  Once  validated  with  model  scale  data, 
computational  ship  hydrodynamics  produces  the  crucial 
damping  information  required  to  predict  full  scale  motions. 

Computational  ship  hydrodynamics  provides  the  tool  to 
explore  the  effects  of  nonconventional  geometry  on  ship  per¬ 
formance.  For  example,  catamaran  and  trimaran  configura¬ 
tions  applied  to  modified  monohulls  for  wave  cancellation, 
and  radical  bow  shaping  for  wave-piercing,  including  sweep 
down  rather  than  sweep  up,  are  approachable  with  computa¬ 
tional  ship  hydrodynamics.  Computational  ship  hydrodynam¬ 
ics  is  already  being  used  to  evaluate  modifications  to  the  bow 
sonar  dome  for  the  DDG-51  as  an  “after  the  fact”  demonstra¬ 
tion.  The  original  modifications  were  made  based  solely  on 
towing  tank  tests,  and  were  limited  by  resources  to  a  selected 
set  of  geometry.  The  goal  with  the  demonstration  is  to  examine 
the  accuracy  of  a  numerical  process,  and  perhaps  to  result  in 
a  better  design  for  the  modification. 

Computational  ship  hydrodynamics  can  be  employed  to 
develop  control  systems,  and  in  particular  by  providing  force  and 
moment  loads  for  input  to  control  system  design  and  optimization 
especially  for  cases  of  extreme  maneuvers.  It  is  very  difficult  to 
obtain  required  towing  tank  (or  more  specifically,  maneuvering 
basin)  model  data  for  firee-running  models.  A  particular  problem 
to  anticipate  by  design  or  operation  extreme  maneuvers  for 
submarines,  or  other  appended  bodies  of  revolution.  Typicalliy 


predictions  are  made  using  system  identification  techniques, 
and  validation  of  the  design  is  attempted  by  measuring  trajec¬ 
tories  for  free-running  models.  Computational  ship  hydrody¬ 
namics  permits  the  scaling  of  the  results  to  full  scale,  and  is  a 
technique  for  analyzing  the  flow  round  the  body  to  determine 
the  causes  of  unfavorable  loadings.  Alternative  physical  meas¬ 
urements  of  the  details  of  the  flow  are  very  expensive  in  both 
money  and  time. 

Another  example  concerns  the  design  of  propellers.  Cur¬ 
rently  they  are  designed  using  a  variety  of  specialized  and 
nonintegrated  computational  methods,  developed  piecemeal 
over  years  of  development.  The  movement  is  toward  flow 
solvers  that  capture  all  the  features,  but  compute  only  what  is 
needed.  For  example,  gridding  can  focus  on  regions  of  the  flow 
of  interest  with  large  density  grids,  with  much  smaller  density 
of  gridding  elsewhere,  to  capture  the  major  effects  on  the 
region  of  interest. 

A  significant  step  towards  a  numerical  towing  tank  will 
soon  take  place  when  “simultaneous”  towing  tank  measur- 
ments  and  RANS  computations  will  be  made  with  the  self-pro¬ 
pelled  DTMB  Model  5415  Surface  Combatant.  The 
computations  and  measurements  will  leap-frog  to  guide  the 
entire  investigation  toward  an  archival  database  of  informa¬ 
tion  on  propeller/hull  interactions,  flow  separations,  wave 
breaking,  and  transom  stem  flows  in  support  of  advanced 
surface  ship  designs. 

Computational  Ship 
Hydrodynamics  Today 

In  the  selected  examples  that  follow,  it  is  shown  that  the 
complicated  three-dimensional  flows  around  various  compo¬ 
nents  of  ships  are  represented  by  computational  ship  hydrody¬ 
namics.  All  of  the  cases  discussed  below  consider  the  turbulent 
flow  around  the  ship  by  including  the  Reynolds  stresses.  The 
flow  severs  are  generally  know  as  “fully  nonlinear  free-sur- 
face  RANS  codes”  or  more  simply  as  “RANS  codes”  in  the 
case  that  that  the  free  surface  is  not  a  factor.  All  of  the 
computational  ship  hydrodynamics  computer  codes  include 
time  in  the  solutions.  In  the  case  of  steady  motions,  evolving 
time  is  used  to  arrive  at  the  steady  motion  and  may  not  be  an 
accurate  representation  of  real  fluid  motion  as  a  trade-off  for 
computational  efficiency.  Although  sophisticated  computa¬ 
tional  methods  that  neglect  turbulence  are  used  for  computa¬ 
tional  fluid  dynamics,  and  will  always  have  a  proper  place  in 
design  methodology  because  of  the  reduced  computational 
resources  required,  the  proper  prediction  of  the  flow  by  com¬ 
putational  ship  hydrodynamics  for  final  design  or  for  analysis 
must  include  turbulence. 

Recent  examples  of  computations  for  surface  ships  in¬ 
clude  the  prediction  of  the  nonlinear  free-surface  turbulent 
flow  around  a  modem  naval  combatant  by  Stem,  et  al.  (1995) 
and  around  a  sailing  yacht  by  Farmer,  et  al.  (1995).  The 
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computation  around  the  surface  combatant  features  the  effects 
of  the  bulbous  bow  and  the  transom  stem.  A  sample  result  from 
the  computation  is  shown  in  Figure  4.  As  with  of  all  the 
advanced  computational  methods,  the  discretization  of  the 
entire  flow  field  permits  the  full  interaction  of  the  low  speed 
boundary  layer  next  the  to  the  hull  and  the  wavemaking  on  the 
free  surface.  The  computations  predict  a  “double  rooster  tail” 
associated  with  the  two  peaks  in  the  free  surface  just  aft  of  the 
stem.  A  comparison  with  measurements  obtained  in  a  towing 
tank  test  shows  that  the  twin  peaks  are  present  in  the  real  flow. 
In  the  second  example  the  computations  for  the  sailing  yacht 
include  the  effects  of  yaw  angle,  in  which  the  flow  around  the 
submerged  keel,  acting  as  a  wing,  generates  a  lifting  force  and 
a  downstream  trailing  vortex  that  interacts  with  the  radder. 
Simulation  and  modeling  of  these  types  of  interactions  aie 
essential  for  the  prediction  of  the  inflow  to  propulsors  and 
mdders,  and  for  the  appendage  interactions  on  maneuvering 
submarines,  especially  near  to  or  broached  at  the  free  surface. 
The  numerical  procedure  that  produced  these  computations  is 
expected  to  be  applied  to  the  America’s  Cup  designs  for  1999 
in  what  promises  to  be  an  intense  effort  by  the  United  States 
to  reclaim  the  trophy. 

Computations  of  the  flow  interaction  between  the  ship 
hull  and  the  operating  propeller  have  been  made  by  Stem,  et 
al.  (1994)  for  a  full  domain  solution  and  by  Weems,  et  al. 
(1994)  for  a  partial-domain  solution  (in  which  the  flow  away 
from  the  hull  is  approximated  by  inviscid  flow)  for  a  ship 
model  used  by  the  International  Towing  Tank  Conference  for 
comparison  of  measurement  and  computation  methods  (Series 
60,  Cb  =  0.6).  In  both  cases  the  propeller  flow  is  computed 
using  potential  flow  methods  that  do  not  fully  account  for  the 
viscous,  turbulent  flow  but  do  predict  the  lifting  force  on  the 
propeller  blades.  This  force  is  embedded  in  the  RANS  com¬ 
putations  as  a  body  force,  distributed  azimuthally  and  radially 
as  appropriate  in  the  propeller  disc  region.  The  result  is  that 
the  interactive  effect  of  the  propeller  and  the  hull  is  captured 
in  the  computations.  This  is  an  important  effect  required  for 
propeller  design;  it  arises  from  the  suction  effect  of  the  pro¬ 
peller  on  the  hull,  which  changes  the  momentum  distribution 
of  the  hull  boundary  layer  upstream  from  the  immediate  region 
of  the  propeller.  The  change  in  the  upstream  boundary  layer 
momentum  then  appears  as  a  change  in  the  effective  open- 
water  inflow  to  the  propeller,  leading  to  a  change  in  propeller 
loading.  Propeller  design  depends  on  the  proper  accounting 
for  this  effective  wake,  which  is  not  measurable.  Comparisons 
of  computed  and  measured  velocity  profiles  for  the  flow 
around  the  propeller  demonstrate  that  the  flow  is  qualitatively 
captured  by  the  computations,  but  that  detailed  resolution  of 
the  flow  is  missing.  The  results  willbe  improved  with  better 
turbulence  models  and  higher  resolution  computational  meth¬ 
ods,  and  these  are  active  research  areas. 

The  marine  propeller  is  a  complex  geometry  with  a  com¬ 
plicated  flow.  Specialized  computational  methods  have  been 


Figure  4. 


Results  of  Computation  of  the  Free-Surface  Deformation  of  the 
How  Around  the  FF1 052  at  Froude  Number-  0.3  and  Reynolds 
Number  4, 000, 000. 


applied  to  the  propeller  and  its  interactions  with  the  flow 
around  the  ship.  Clearly  modem  propellers  with  large  three- 
dimensional  curvatures  to  the  blades  and  possibly  with  a 
surrounding  duct  are  complex  geometries.  Specification  of 
these  geometries,  and  the  maintenance  of  the  surface  integrity 
in  gridding  of  the  surface  is  a  specialized  subject.  The  details 
of  the  propeller,  such  as  the  intersection  of  the  blade  root  with 
the  hub  and  the  clearance  between  the  blade  and  the  duct,  must 
be  precisely  specified  in  cases  where  the  vortical  flows  in  such 
regions  can  lead  to  nuisance  cavitation.  The  complex  geome¬ 
try  produces  a  complex  flow  which  is  made  more  difficult  to 
predict  by  the  interactions  between  the  propeller  flow  and  the 
flow  around  the  hull  as  described  above.  The  need  for  even 
higher  performance  propellers  has  led  to  the  need  to  predict 
the  hull  and  propeller  flow  simultaneously.  Kerwin,  et  al. 
(1994)  have  developed  a  design  method  for  complex  propul¬ 
sors  which  includes  a  rational  B-spline  method  for  geometry 
specification  as  well  as  a  blend  of  RANS  computations  for  the 
axisymmetric  inflow  and  a  lifting  surface  method  for  the  blade 
loading.  Figure  5,  shows  a  geometry  characteristic  of  modem 
propulsors.  There  are  very  significant  flow  interactions  be¬ 
tween  the  components,  including  the  blade-to-blade  and 
blade-row-to-row,  and  the  upstream  boundary  layer  and  the 
flow  through  the  propulsor.  Such  designs  are  further  compli¬ 
cated  by  the  full  integration  of  the  propulsor  and  the  stem 
appendage  flows.  With  this  type  of  integration  the  stem  and 
appendages  and  the  propulsor,  considered  separately,  will  not 
meet  requirements.  But  the  combination,  working  together, 
does  meet  the  requirements.  To  compute  the  flow  for  these 
types  of  arrangements  it  is  necessary  to  consider  the  full 
turbulent  flow  through  the  system. 

Significant  and  substantial  progress  has  recently  been 
made  in  the  computation  of  the  unsteady  turbulent  flow  around 
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Figure  5. 

Drawing  of  the  Complex  Geometry  Characteristic  of  Modern 
Turbomachinery-Like  Propulsors. 


a  propeller.  (Stem,  et  al.,  1994)  describe  a  computational 
method  for  steady  flow  around  a  marine  propeller  with  a  large 
amount  of  curvature  to  the  blades.  Such  a  highly  three-dimen¬ 
sional  geometry  produces  complicated  three-dimensional 
flow  separations,  which  contribute  to  blade  passage  vortices 
and  the  wake  of  the  propeller.  The  computational  method  was 
validated  with  experimental  data  and  found  to  be  at  least  as 
accurate  as  the  existing  inviscid  computational  methods.  The 
added  capability  of  this  RANS  method  is  the  prediction  of  the 
turbulent  viscous  wake  without  approximation,  except  for  the 
turbulence  model.  It  is  expected  that  such  methods,  free  of  the 
artificial  models  incorporated  in  inviscid  models,  will  accu¬ 
rately  treat  complicated  geometry  such  as  that  shown  in  Figure 
5.  A  companion  paper  by  Chen,  et  al.  (1994)  describes  a 
successful  attempt  to  compute  the  unsteady  flow  around  a 
lifting  surface.  This  method  is  being  extended  to  the  complete 
flow  around  a  propeller  operating  in  a  spatially  nonuniform 
inflow. 

The  computation  of  flows  around  submarines  generally 
excludes  the  need  to  include  the  free  surface,  but  necessarily 
includes  increased  emphasis  on  predicting  vortical  flows. 
These  vortical  flows  arise  as  the  lifting  surfaces  used  for 
control  generate  trailing  vortices  both  at  their  root  and  tip  ends 
as  the  consequence  of  lift,  and  at  their  root  in  addition  as  an 
interaction  with  the  hull  boundary  layer,  and  along  the  hull  if 
it  is  at  incidence  to  the  flow.  These  vortical  flows  interact  with 
downstream  surfaces,  such  as  other  appendages  and  the  pro¬ 
peller.  The  computational  burden  is  to  accurately  predict  the 
evolution  of  the  vortical  flow,  which  is  characterized  by  large 
velocity  gradients  and  the  suppression  of  turbulence  in  the 


core.  Computational  methods  must  not  be  overly  dissipative 
or  the  vorticity  will  artificially  dissipate;  hence  the  computer 
codes  are  very  sensitive  to  moderate  changes  in  input. 

RANS  computational  methods  have  been  successfully 
used  to  model  the  separated  flow  around  a  body  of  revolution 
at  an  angle  of  incidence.  This  flow  is  complicated  by  the 
three-dimensional  flow  separations  that  occur  along  the  sides 
of  the  body.  Such  separations  are  difficult  to  model  because 
the  flow  detaches  from  the  smooth  surface  of  the  body  rather 
than  from  a  sharp  edge  defining  the  point  of  detachment.  Sung 
&  Huang  ( 1 994)  have  produced  validated  computations  for  the 
forces  and  moments  on  the  body  at  incidence  for  large  model- 
scale  Reynolds  numbers.  In  order  to  achieve  the  successful 
computations,  it  was  necessary  to  customize  the  turbulence 
model  for  the  appearance  of  the  longitudinal  vortices  gener¬ 
ated  by  the  flow  separations.  Submarines  have  appendages, 
which  must  be  included  in  the  full  computation  of  the  flow 
around  the  submarine.  A  chief  difficulty  is  the  prediction  of 
the  flow  around  the  intersection  of  the  appendage  with  the 
body.  Such  intersections  produce  a  roll-up  of  the  boundary 
layer  fluid  along  the  body  to  produce  horseshoe  vortices 
wrapping  around  the  leading  edge  of  the  appendage.  The 
consequence  is  that  the  boundary  layer  is  inverted  downstream 
from  the  appendage.  This  feature  must  be  captured  by  high 
spatial  resolution  computations.  An  efficient  solver  for  this 
flow  has  been  validated  with  experimental  data  (Sung  & 
Huang,  1994).  The  solver  utilized  multiblock,  multi  grid,  local 
refinement,  preconditioning,  and  adaptive  artificial  dissipa¬ 
tion  models  to  achieve  90  percent  savings  in  both  memory  and 
cpu  time.  Such  efficiencies  are  required  to  permit  practical 
computational  ship  hydrodynamics  applications. 

The  research  group  at  the  Applied  Research  Labora¬ 
tory/Pennsylvania  State  University  has  been  using  unsteady 
RANS  flow  solvers  for  unsteady  motion  to  determine  the 
forces  and  moments  on  a  submarine  (see  Figure  6).  With  this 
information,  predictions  of  the  trajectory  of  the  submarine  are 
possible.  This  work  is  ongoing,  and  it  is  expected  that  the 
computations  will  be  used  to  design  a  validation  test  in  which 
the  computed  and  measured  trajectories  will  be  compared  in 
the  near  future.  Current  research  with  surface  ships  is  also 
emphasizing  time-dependent  solutions.  Particular  interest  is 
focused  on  turns  in  shallow  water,  and  the  effects  of  channel 
crossings. 

Surface  ships  are  characterized  by  motions  induced  by 
interactions  with  waves  as  well  as  by  maneuvers  of  the  vehicle 
itself.  Recent  efforts  to  predict  ship  motions  include  the  effects 
of  curved  hulls  in  the  transverse  cross-section  and  the  transom 
stem.  These  geometry  features  have  important  effects  on  the 
motions  and  loads  on  ships  because  they  are  major  factors  in 
the  hydrostatic  loading  and  the  viscous  flow.  The  state-of 
-the-art  for  ship  motion  computations  is  inviscid  panel  meth¬ 
ods.  These  methods  rely  on  various  boundary  integral  formu¬ 
lations  to  represent  the  flow  around  the  ship.  Viscous  effects 
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Figure  6. 

Computed  Flow  Trajectories  for  a  Submarine-Like  Body  Un¬ 
dergoing  Sudden  Pitch  (Whitfield,  1995b). 


are  accounted  for  at  the  stem  by  enforcing  a  Kutta  condition 
so  that  the  flow  smoothly  exits  from  the  hull.  Sclavounos, 
Nakos  &  Huang  (1994)  have  developed  a  robust  computa¬ 
tional  method  (SWAN)  that  treats  the  unsteady  frequency 
domain  ship  motion  problem  as  a  linear  perturbation  of  the 
nonlinear  solution.  SWAN  accounts  for  tiiree-dimensional 
flow  effects  for  hulls  with  flare  and  transom  stems,  and  pre¬ 
liminary  results  show  very  good  agreement  with  experimental 
data  for  pitch  and  heave  motions.  The  inviscid  methods  are 
more  suited  for  pitch  and  heave  than  for  roll  and  yaw,  which 
are  dependent  on  viscous  effects  along  the  entire  hull  and 
require  a  full  treatment  of  viscosity  and  turbulence.  Neverthe¬ 
less,  the  inviscid  SWAN  method  is  currently  under  validation 
studies  as  a  significant  improvement  to  the  current  two-dimen¬ 
sional  methods  used  for  traditional  ship  motion  design.  A 
significant  research  effort  is  underway  to  develop  and  apply 
unsteady  RANS  computations  to  the  rollling  motion  of  a 
surface  ship  in  a  seaway.  The  damping  of  roll  motion  is 
strongly  dependent  on  the  viscous  turbulent  flow  and  cannot 
be  represented  by  inviscid  methods. 

Numerical  Simulations  and 
Surrogate  Models 

Computational  ship  hydrodynamics  is  based  on  solutions 
to  the  unsteady  RANS  equations.  It  is  not  possible  with  those 
equations  to  explicitly  compute  all  the  details  of  the  flow.  The 
equations  themselves  are  averaged  over  a  short  characteristic 
time,  a  process  that  inherently  loses  detail.  A  more  practical 
problem  arises  in  consideration  of  the  computational  resources 
in  both  time  and  hardware  that  would  be  required  to  compute 
the  details  for  short  characteristic  times  as  well  as  the  longer 
ship  motion  time  scales.  In  practice,  models  are  used  to  repro¬ 
duce  the  effects  of  the  small  details  of  the  flow  on  the 
macroscale  features.  These  models  can  be  obtained  through 
interpretation  of  physical  measurements  or  of  direct  numerical 


simulations  of  the  flow.  Numerical  simulations  are  often  pre¬ 
ferred  because  spatial  and  temporal  resolution  can  be  greater 
than  in  a  physical  measurement,  and  because  many  variables 
can  be  computed  which  are  not  easily  measured  (such  as  the 
pressure-velocity  correlation).  As  much  as  computational 
ship  hydrodynamics  is  not  complete  without  the  performance 
validations  at  the  end,  it  is  also  not  complete  without  the 
simulations  and  modeling  required  at  the  beginning  to  develop 
the  equations  to  be  solved.  There  are  two  examples  considered 
here,  the  simulation  of  the  disconnection  of  the  vortex  ring  as 
a  turbulence  mechanism  and  the  modeling  of  the  turbulent 
wake  of  the  spilling  breaking  wave  as  a  RANS  surrogate 
model. 

Vortex  Disconnection 

Free-surface  turbulence  models  are  needed  to  represent 
the  evolution  of  the  turbulent  flow  near  the  free  surface.  In  this 
case  the  boundary  is  shear- free,  and  is  different  than  the  no-slip 
boundary  characteristic  of  the  turbulent  flow  along  the  ship 
hull.  A  fundamental  mechanism  for  free-surface  turbulence  is 
the  vortex  disconnection  event  during  which  vorticity  appears 
(or  disappears)  on  the  free-surface  boundary  of  the  fluid.  This 
process  is  analogous  to  the  hairpin  vortex  mechanism  under¬ 
lying  solid  body  boundary  layer  mechanics.  Observations  of 
the  event  were  first  made  for  a  vortex  ring  interacting  with  the 
free  surface  (Bernal  &  Kwon,  1989);  a  mechanism  was  pos¬ 
tulated  by  Rood  (1994)  using  a  rigorous  analysis  of  the  mathe¬ 
matical  physics;  and  laboratory  measurements  confirmed  the 
presence  of  the  mechanism  (Gharib,  1994).  The  details  of  the 
three-dimensional  flow  as  it  evolves  in  time  were  investigated 
with  numerical  experiments,  which  proved  to  be  the  key  to 
obtaining  the  needed  quantification  of  the  various  dynamics 
underlying  the  event  (Yue,  1995).  Figure  7  is  a  computer-gen¬ 
erated  picture  of  one  stage  of  the  event.  Because  the  computer 
computations  have  captured  the  quantitative  information  for 
all  the  variables,  the  computations  can  be  used  to  obtain  details 
important  for  assessing  the  role  of  various  terms  in  turbulence 
models. 

Breaking  Spilling  Wave 

An  example  flow  that  must  be  modeled  for  computational 
ship  hydrodynamics  is  the  wake  of  the  spilling  breaking  wave, 
such  as  the  bow  wave.  A  scenario  for  the  model  for  the  wake 
of  the  spilling  breaking  wave  consists  of  several  embedded 
numerical  computations.  Initially  the  RANS  computation  of 
the  free-surface  deformation  indicates  that  the  wave  is  build¬ 
ing  to  a  height  indicative  of  incipient  wave  breaking  (Paterson, 
Tahara,  &  Stem,  1995).  A  numerical  simulation  using  exact 
physics  is  performed  to  determine  the  onset  of  wave  breaking 
(Dommermuth  &  Mui,  1995).  The  breaking  wave  is  then 
modeled  as  a  mixing  layer  with  vorticity  flux  emanating  from 
the  toe  of  the  wave  (Cointe  &  Tulin,  1994).  The  value  for  the 
vorticity  flux  is  obtained  from  the  local  velocity  upstream  from 
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Figuro  7. 

Numerical  Simulation  of  Vortex  Loop  Disconnecting  to  the  Free 
Surface  (Successive  Time  from  Top  to  Bottom). 


the  wave  as  determined  from  the  RANS  computation  of  the 
three-dimensional  nonbreaking  wave  (Rood,  1995).  A  large 
eddy  simulation  is  employed  to  compute  the  roll  up  and 
production  of  the  turbulent  features  of  the  mixing  layer,  which 
are  then  characterized  as  a  turbulence  spectrum.  This  turbu¬ 
lence  spectrum  is  used  to  compute  the  Reynolds  stresses, 
which  are  then  inserted  as  a  turbulence  model  in  the  RANS 
computations  immediately  downstream  from  the  breaking 


wave.  The  RANS  computations  then  predict  the  evolution  of 
the  turbulent  wake  as  it  moves  downstream  from  the  wave. 

The  turbulence  model  for  the  breaking  wave  is  a  surrogate 
model  that  is  obtained  from  embedded  simulations  within  the 
RANS  computation.  The  surrogate  model  links  the  informa¬ 
tion  from  the  detailed  physics  obtained  with  the  embedded 
simulation  to  the  macroscale  flow  description  (or  effects) 
required  to  advance  the  hydrodynamic  system  computations. 
The  embedded  simulations  are  computationally  burdensome, 
and  in  this  scenario  they  are  made  only  occasionally  in  time 
and  locally  in  space  as  needed  to  obtain  the  coefficients 
required  for  the  turbulence  model  for  breaking  wave. 

The  Future 

The  emergence  of  computational  ship  hydrodynamics  is 
a  beginning.  It  is  the  beginning  of  a  major  change  in  the  way 
ships  and  submarines  are  designed,  analyzed,  and  operated. 
The  issue  is  how  to  manage  and  capitalize  on  this  revolution 
to  make  it  affordable.  Intelligent  and  informed  leadership  is 
required  to  maintain  a  steady  helm  and  an  even  keel  along  the 
path  to  simulation-based  design,  the  ultimate  goal. 

Computational  ship  hydrodynamics  is  not  inexpensive.  A 
typical  RANS  computation  for  the  flow  around  a  surface 
combatant  currently  consumes  some  70  CPU  hours  on  a  su¬ 
percomputer.  To  be  a  useful  engineering  tool,  computational 
ship  hydrodynamics  must  suffer  the  processes  of  accreditation 
and  certification.  Naval  hydrodynamicists  must  assume  that 
computer  hardware  capability  and  algorithm  development  will 
pace  each  other  to  produce  orders  of  magnitude  more  capabil¬ 
ity. 

Is  computational  ship  hydrodynamics  more  (or  less)  ex¬ 
pensive  than  the  cost  for  towing  tank  tests?  This  is  like  com¬ 
paring  apples  and  oranges;  computations  and  measurements 
must  be  viewed  as  complementary  tools  with  an  overlap  in 
which  validation  is  obtained.  In  any  case,  towing  tank  meas¬ 
urements  themselves  will  become  more  costly  and  sophisti¬ 
cated  as  such  innovations  as  digital  particle  image  velocimetry 
are  employed  to  measure  the  time  dependent  flow  of  large 
eddies  around  ship  hulls  to  develop  higher  performance  ships. 
Computations  cannot  practically  capture  all  the  details  of  the 
flow,  which  contribute  to  the  macroscale  features  such  as  the 
breaking  wave  and  air  entrainment.  Measurements  cannot 
capture  all  the  parameterization  available  to  computations, 
such  as  ranges  of  angle  of  incidence.  Once  the  computation  is 
established  (converged  and  stable,  and  validated  and  cali¬ 
brated),  the  geometry  can  be  modified  with  ease  to  produce  a 
range  of  information. 

At  the  turn  of  the  century,  naval  hydrodynamicists  expect 
computational  ship  hydrodynamics  to  predict  the  hydrody¬ 
namic  flow  system,  consisting  of  the  hull,  appendages,  propel¬ 
ler  for  arbitrary  ship  orientation  in  a  single  frequency  seaway. 
The  computations  will  predict  maneuvering  and  transient  shal- 
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low  water  operations.  Such  a  code  will  be  a  numerical  un¬ 
steady  Reynolds- Averaged  Navier  Stokes  prediction  method. 
The  computer  code  will  be  validated  but  may  require  calibra¬ 
tion  against  towing  tank  test,  and  will  be  used  to  extend  the 
range  of  parameters  and  to  scale  results  to  full  scale. 
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Turbulence:  A  Roadblock 
to  Computational  Fluid 
Dynamics 

L.  Patrick  Purtell 
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Abstract 

Continuing,  indeed  accelerated,  advances  in  computa¬ 
tional  fluid  dynamics  (CFD)  are  required  to  meet  the  chal¬ 
lenges  of  simulation-based  analysis  and  design  of  naval 
hydrodynamics  systems.  The  great  progress  achieved  in  the 
power  of  computers  and  in  numerics  applicable  to  the  equa¬ 
tions  of  fluid  mechanics  exposes  the  critical  need  for  adequate 
representations  of  turbulence.  An  existing  extensive  hierarchy 
of  turbulence  models  must  be  evaluated,  and  new,  more  so¬ 
phisticated  models  need  to  be  developed  for  the  truly  complex 
flow  fields  of  practical  interest. 

Introduction 

Richard  Feynman  is  said  to  have  remarked,  “Turbulence 
is  the  last  great  unsolved  problem  of  classical  physics.”  De¬ 
spite  decades  of  extensive,  dedicated,  innovative,  even  ingen¬ 
ious  attempts  to  derive  a  “theory  of  turbulence”  which  can 
predict  its  statistical  characteristics  (eg.,  average  fluctuation 
levels),  turbulence  as  a  whole  remains  beyond  our  analytical 
capability.  Nevertheless  there  is  a  very  practical  need  to  pre¬ 
dict  turbulent  flowfields  in  remarkably  diverse  conditions: 
around  ships  and  aircraft,  in  pipes  and  ducts,  in  chemical 
mixers  and  combustors,  and  in  the  atmosphere  and  ocean. 

This  need  to  provide  solutions  to  pressing  technical  prob¬ 
lems  led  to  the  development  over  many  years  of  numerous 
simple  correlations  and  approximate  ad  hoc  methods  in  fluid 
mechanics.  However,  the  explosive  growth  in  computer  power 


and  associated  numerical  techniques  in  the  last  few  decades 
has  correspondingly  reduced  the  need  for  such  methods.  The 
ability  to  calculate  discretized  versions  of  the  equations  of 
fluid  mechanics  has  evolved  rapidly  and  has  been  limited 
primarily  by  computer  power.  The  representation  of  turbu¬ 
lence  in  these  codes,  however,  has  not  kept  pace,  though 
significant  advances  have  been  made.  As  computer  power 
increases  and  numerical  methods  are  improved,  the  limitations 
to  accuracy  caused  by  inadequate  turbulence  models  become 
more  and  more  apparent. 

A  brief  summary  of  the  fundamentals  underlying  contem¬ 
porary  turbulence  models  is  in  order  here.  The  Navier-Stokes 
Equations  (NSE),  the  fundamental  equations  of  viscous  flow 
relate  changes  in  momentum  to  forces 

3m,-  3m,-  1  3a', y 

where  U;  is  the  velocity  and  Oy  is  the  stress  tensor.  Two 
characteristics  of  the  NSE  in  particular  should  be  noted:  first, 
the  nonlinearity  of  the  second,  “convective”,  term  (the 
source  of  “richness”  for  the  physicist  and  of  “problems”  for 
the  engineer!)  and  second,  the  requirement  that  at  a  solid 
boundary  the  fluid  velocity  equals  the  velocity  of  the  bound¬ 
ary  (eg.,  zero  for  a  stationary  boundary).  For  flows  of  low 
momentum  and  high  viscosity  (ie.,  flows  of  low  Reynolds 
number,  the  ratio  of  momentum  to  viscous  resistance),  the 
flow  will  be  smooth  and  generally  predictable,  even  in  un¬ 
steady  conditions.  Numerous  analytical  solutions  to  the  NSE 
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have  been  derived  for  such  conditions,  and  today  computa¬ 
tional  techniques  render  such  flows  effectively  solved.  How¬ 
ever,  such  flows  are  rare  in  practice,  and  at  realistic  Reynolds 
numbers,  the  flow  is  unstable  and  becomes  turbulent 
through  a  variety  of  dynamical  processes  (the  study  ofwhich, 
'‘stability  and  transition”,  is  a  field  of  its  own). 

Turbulence  has  numerous  intriguing  characteristics 
which  would  require  volumes  to  describe  in  detail.  For  present 
purposes  it  will  suffice  to  summarize  only  a  few  of  them.  First, 
turbulence  has  not  been  found  to  disobey  the  NSE  at  any  point 
or  time  in  the  flow.  This  is  crucial,  of  course,  since  otherwise 
we  would  need  to  rederive  the  fundamental  equations  of 
motion.  Second,  turbulence  involves  fluctuations  over  a  wide 
range  of  length  and  time  scales  ^  a  range  which  increases  with 
the  Reynolds  number.  The  largest  fluctuation,  or  “eddy”,  is  of 
the  order  of  the  size  of  the  flow,  such  as  the  width  of  a  wake 
or  the  diameter  of  a  pipe.  The  smallest  scale  is  limited  by 
viscous  dissipation,  but  it  is  typically  orders  of  magnitude 
smaller  than  the  largest  scales. 

Third,  the  eddies  are  strongly  affected  by  the  average 
motion  and  likewise  affect  the  average  motion.  They  are  not 
merely  “swept  along”  with  the  flow.  It  should  be  noted, 
though,  that  turbulent  flow  is,  on  average,  “similar  to”  laminar 
flow;  it  is  not  radically  different  in  overall  character  (eg.,  a 
wake  is  still  a  wake).  This  is  critical  for  computational  ap¬ 
proaches.  Fourth,  and  finally,  turbulence  is  important  It  can¬ 
not  be  ignored,  and  it  cannot  often  be  replaced  by  simple 
approximations. 


Representations  of  Turbulence 

Observation  of  eddies  and  thus  a  perception  of  turbulence 
long  predates  the  modem  development  of  fluid  mechanics,  so 
it  is  no  surprise  that  attempts  to  account  for  turbulence  analyti¬ 
cally  nearly  coexisted  with  the  development  of  the  NSE  them¬ 
selves.  The  gross  “similarity”  of  turbulent  flows  to  laminar 
flows  mentioned  before  soon  prompted  the  concept  of  an 
“eddy  viscosity”  analogous  to  the  molecular  viscosity.  Thus 
the  expression  for  the  shear  stress  in  unidirectional  laminar 
flow,  o' which  is  proportional  to  the  rate  of  strain  ( v  is  the 
kinematic  viscosity  and  p  is  the  density): 


0'l2  =  pv 


dxn 


is  replaced  by 


572-P(v  +  v,)-j^ 

where  the  overbar  indicates  a  time-average  (or  ensemble 
average  for  unsteady  flows),  and  an  “eddy  viscosity”,  u  j,  has 
been  added  to  the  molecular  kinematic  viscosity.  Dimension- 
ally,  the  eddy  viscosity  is  a  characteristic  length  (squared) 


overa  characteristic  time.  A  straightforward  use  ofa  constant 
eddy  viscosity  would  imply  that  the  turbulence  is  a  property 
of  the  fluid,  whereas  it  is  a  strong  property  of  the  flow,  varying 
widely  throughout  a  typical  flowfield.  Simple  expressions  for 
the  eddy  viscosity  are  at  the  foundation  of  many  turbulence 
models,  but  work  only  in  flows  for  which  they  are  tightly 
calibrated  and  in  conditions  where  the  turbulence  is  in 
approximate  equilibrium  with  the  mean  flow,  eg.,  in  far 
wakes  or  jets.  It  has  been  pointed  out^  that  “...  eddy  viscosity 
concepts  underlie  the  entire  range  of  models  for  turbu¬ 
lence...”  from  representing  all  of  the  turbulence  (simple 
models)  to  representing  specific  quantities  (such  as  terms  in 
higher  order  transport  models).  The  latter  are  really  gener¬ 
alized  eddy  viscosities  which  may  vary  in  space  and  time. 

There  are  several  broad  categories  of  theoretical  ap¬ 
proaches  to  turbulence,  eg.,  statistical  theories  utilizing  the 
extensive  body  of  mathematics  developed  for  statistical  me¬ 
chanics,  wave  (Fourier)  descriptions  concerned  in  particular 
with  energy  flow  among  scales  of  motion,  and  more  recently 
dynamical  systems  theory  (chaos).  These  will  not  be  discussed 
here  where  the  focus  is  on  the  more  utilitarian  aspects  of 
modeling  and  its  use  in  computational  fluid  dynamics  (CFD). 
However,  advances  in  fundamentals  of  modeling  are  likely  to 
depend  more  and  more  on  sound  theoretical  concepts.  Rational 
approaches  utilizing  renormalization  group  theory  or  rapid 
distortion  theory  are  more  frequently  found  in  the  literature 
today  underlying  important  advances  in  sophisticated  models. 

Osborne  Reynolds^  introduced  a  concept  which  today 
forms  the  basis  of  the  representation  of  turbulence  in  most  of 
computational  fluid  dynamics  (CFD).  He  decomposed  the 
velocity  into  a  time-  (or  ensemble-)  averaged  part,  Uj,  and  a 
fluctuating  part,  Uj. 


u'i  =  Ui  +  Ui 

He  then  introduced  this  decomposition  into  the  equations 
of  motion,  averaged  the  equations  over  time  (or  ensembles), 
and  produced  what  today  are  named  the  Reynolds-Averaged 
Navier-Stokes  (RANS)  equations.  The  nonlinearity  of  the 
convective  terms  produces  an  average  (or  correlation)  of  the 
product  of  the  fluctuating  components  termed  the  “Reynolds 
stress”,  ufuj 


du\ 

u':-=r^-^U: 


dX: 


dUi  du^Uj 
^  dXj  dxj 


where  u^Uj  is  frequently  denoted  by  t^.  The  eddy  viscosity 

concept  then  has  t-  set  proportional  to  the  average  rate-of- 
strain. 

Simple  expressions  of  this  concept  wherein  the  eddy 
viscosity  is  determined  by  algebraic  expressions  are  termed 
“zero-equation”  models  since  no  auxiliary  differential  equa¬ 
tion  must  be  solved.  The  more  developed  examples  of  these, 
eg.  Cebeci-Smith  and  Baldwin-Lomax  models,  are  popular  for 
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boundary  layers  because  of  their  simplicity  and  low  demands 
on  computer  resources.  They  are,  however,  inaccurate  for 
complex  flows  and  require  “calibration”  for  the  particular  type 
of  flow  encountered. 

Instead  of  employing  algebraic  expressions  for  the  char¬ 
acteristic  length  and  velocity  (and  thus  the  eddy  viscosity), 
“moments”  of  the  NSE  may  be  derived.  The  NSE  may  be 
considered  as  “transport  equations”  for  the  momentum.  If  they 
are  multiplied  by  velocity  and  then  averaged,  a  transport 
equation  for  kinetic  energy  may  be  derived,  and  in  a  similar 
manner  transport  equations  for  other  quantities  may  be  devel¬ 
oped. 

This  approach  demands  some  further  discussion.  The 
higher  moment  equations  are  derived  completely  from  the 
NSE  and  thus  contain  no  additional  information  about  the 
flowfield!  It  is  merely  hoped  that  by  modeling  terms  in  the 
higher  order  equations  instead  of  in  the  first  order  we  may 
incorporate  more  complex  physics  in  the  flow.  This  can  indeed 
be  convincingly  argued,  but  proper  models  for  some  of  the 
complicated  higher  order  terms  are  so  difficult  to  construct  and 
evaluate  that  the  cumulative  errors  may  overwhelm  the  advan¬ 
tages  of  more  physics. 

Models  which  utilize  a  transport  equation  for  the  turbulent 
kinetic  energy  alone  are  the  “one-equation”  models  in  which 
the  eddy  viscosity  is  written: 

v,  =  K‘'^lo 

where  is  a  length  characteristic  of  the  turbulent  eddies  and 
K  is  the  turbulent  kinetic  energy.  Since  we  must  still  devise 
a  method  for  specifying  these  models  are  considered  by 
some  to  be  not  much  of  an  improvement  over  the  zero-equa¬ 
tion  models.  (There  are  exceptions  for  specific  cases,  and  they 
are  still  much  simpler  to  incorporate  in  numerical  codes 
than  the  second  moment  classes  of  models.) 

If  a  second  transport  equation  is  written  for  a  quantity 
which  determines  typically  e ,  the  viscous  dissipation  of 
turbulence,  then  we  have  the  “two-equation”  class  of  models. 
These  models  are  considered  the  simplest  “complete”  models 
in  that  both  the  characteristic  length  and  the  characteristic 
velocity  are  determined  from  PDE’s  and  not  ad  hoc  functions. 
Sophisticated  extensions  to  include  anisotropic  eddy  viscosi¬ 
ties  or  other  complexities  are  still  being  investigated. 

An  alternative  to  an  eddy  viscosity  formulation  involves 
writing  a  transport  equation  for  itself.  The  resulting  “second 
moment”  equations,  or  “second  order  closure”  (SOC),  for  the 
Reynolds  stress  components  comprise  the  most  complex  rep¬ 
resentation  of  turbulence  which  can  be  modeled  for  practical 
purposes.  Because  of  the  nonlinearity  of  the  NSE,  however, 
third  order  terms  appear  in  these  equations  which  must  be 
modeled.  Higher  order  equations  would  likewise  contain  even 
higher  order  terms;  this  is  referred  to  as  the  “closure”  problem 
in  turbulence. 


Variations  on  the  standard  RANS  theme  arise  from  time 
to  time  which  have  intriguing  characteristics.  One  such  ap¬ 
proach  which  gets  occasional  mention  is  sometimes  termed 
the  “split-spectrum  method”.'*  This  technique,  and  others  simi¬ 
lar  to  it,  prescribe  multiple  time  or  length  scales  (not  just  the 
single  scales  of  standard  RANS).  Because  of  increased  com¬ 
plexity  and  added  difficulty  of  calibrating  such  models,  they 
have  not  been  extensively  pursued.  However,  the  availability 
of  direct  simulation  data  today  and  the  realization  that  more 
sophisticated  models  will  be  needed  for  complex  flows  may 
bring  a  revival  of  these  and  other  approaches. 

Though  the  focus  here  is  on  turbulence  modeling  as 
related  to  CFD,  an  important  topic  which  has  the  potential  to 
describe  some  problems  not  addressable  by  RANS  deserves 
mention.  Large-eddy  simulation^,  wherein  the  smaller  scales 
of  turbulence  are  modeled  and  the  larger  scales  are  directly 
computed,  is  a  topic  of  intense  development.  It  is  intermediate 
between  direct  numerical  simulation  (DNS)  of  the  NSE  with 
its  tremendous  demand  on  computer  resources  and  RANS 
which  describes  turbulence  only  at  one  point.  One  particular 
need  in  naval  applications  is  the  ability  to  compute  spatial 
correlations  of  pressure  fluctuations  on  surfaces  for  acoustic 
source  predictions.  Models  for  the  small  (subgrid)  scales  are 
developed  utilizing  the  same  theories  underlying  RANS  de¬ 
velopment,  but  the  hope  is  that  these  small  scales  are  much 
more  universal  in  nature  than  turbulence  as  a  whole  and  thus 
capable  of  description  over  a  wide  range  of  flow  conditions. 

This  concludes  a  very  brief  and  coarse-grained  summary 
of  turbulence  modeling  as  pursued  today.  The  models  are 
numerous  and  vary  in  level  of  complexity,  details  of  calibra¬ 
tion,  and  numerical  implementation.  Their  performance  like¬ 
wise  varies  on  a  wide  range  of  flow  configurations  which  have 
been  addressed. 

CFD:  Equations,  Numerics,  and  Grids 

Computational  Fluid  Dynamics  (CFD)  is  the  prediction 
of  fluid  flow  by  numerical  solution  of  equations  of  motion. 
Equations  are  selected,  discretizations  and  numerical  proce¬ 
dures  are  chosen,  and  the  geometry  of  the  flowfield  specified. 
All  of  these  are  important  and  difficult  steps  in  constructing  a 
numerical  model  of  the  flow  under  consideration.  However, 
several  decades  of  intensive  study  have  raised  many  aspects 
of  the  art  to  a  science  so  that  today  CFD  is  routinely  used  in 
the  design  process  of  many  fluid  flow  systems. 

Only  fairly  recently  has  computer  power  become  so  avail¬ 
able  and  inexpensive  that  solution  of  the  full  NSE  have  been 
utilized  for  design  purposes.  For  large  and  complex  systems 
where  viscous  effects  may  be  (or  at  least  are  declared  to  be) 
unimportant  or  accounted  for  in  other  ways,  inviscid  methods 
are  utilized.  These  may  be  derived  in  terms  of  a  velocity 
potential  (“potential”  methods)  or  by  setting  the  viscosity  to 
zero  in  the  NSE  (“Euler”  methods). 
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Appropriate  numerical  techniques  and  corresponding 
grids  are  then  selected  for  the  computations.  The  importance, 
and  difficulty,  of  these  two  steps  should  not  be  underestimated. 
However,  there  have  been  tremendous  advances  in  these  areas, 
both  in  reducing  required  computer  resources  and  in  improv¬ 
ing  accuracy.  Thus  difficulties  in  turbulence  modeling  when 
using  a  Navier-Stokes  code  become  more  apparent  and  more 
important  to  further  improvements. 

Performance 

One  difficulty  in  evaluating  a  turbulence  model  is  the 
separation  of  errors  due  to  the  model  from  errors  due  to  other 
sources  such  as  the  numerics,  the  grid  resolution,  and  perhaps 
the  most  difficult,  the  interaction  of  the  model  with  these 
others.  However,  techniques  for  estimating  numerical  error  by 
analysis  and  for  determining  grid  problems  by  grid  refinement 
studies  can  in  general  isolate  sources  of  problems.  One  other 
problem  which  must  be  considered,  though,  is  that  of  error  in 
experimental  data  used  for  comparison  with  calculations. 
Here,  too,  there  are  techniques  for  error  analysis  if  they  are 
employed  properly. 

Prediction  of  viscous  flow  around  surface  ships  has  re« 
ceived  more  attention  recently  than  in  the  past  and  has  brought 
to  light  the  difficulties  in  the  computation  of  this  highly 
complex,  vortical,  high  Reynolds  number  flow.  One  effort^ 
which  included  the  needed  examination  of  error  sources  con¬ 
cluded,  “Grid  refinement  in  the  near  wake  or  application  of 
new  accurate  discretization  schemes  have ...  little  effect  on  the 
quality  of  the  solution  in  the  wake.”  For  their  particular  turbu¬ 
lence  models,  a  zero-equation  model  and  several  two-equation 
k-G  models,  the  authors  concluded  that  the  errors  were  caused 
by  “...  too  high  levels  of  turbulence  viscosity.”  They  provided 
further  evidence  of  this  by  arbitrarily  decreasing  the  eddy 
viscosity  in  the  critical  hull  vortex  region  and  got  much  better 
agreement  with  experimental  data.  (See  also  the  companion 
article  in  this  volume,  “The  Emergence  of  Computational  Ship 
Hydrodynamics”) 

The  design  of  aircraft  has  long  been  a  driver  for  advances 
in  CFD.  Here,  as  in  any  engineering  environment,  “...different 
accuracy  demands  (in  design  work)  lead  to  different  accuracy 
demands  on  the  design  tools.”^  The  authors  of  that  review 
paper  also  concluded  that  turbulence  models  were  critical  to 
meeting  design  requirements:  “In  general  weakly  interacting 
flows  can  be  predicted  well...  [but]  ...  strongly  interacting 
flows  of  all  kinds...  often  cannot  be  predicted  well.”  In  particu¬ 
lar  they  propose  focusing  turbulence  modeling  on  complex 
flows  exemplified  by  three-dimensional  separation. 

Some  flows  are  more  strongly  influenced  by  turbulence 
than  others.  A  recent  study  of  flow  over  a  wing-body  configu¬ 
ration^  found  that  for  various  low-order  turbulence  models, 
“...the  wing  pressure  distribution  is  predicted  very  well  [at  zero 
lift]...”.  However,  the  drag  coefficient  was  overpredicted  by 


Figure  1. 

Computed  and  experimental  pressures  atSCP  angle  of  attack. 


nearly  30%.  Also,  the  lift  coefficient  was  predicted  only  to 
within  5%  for  a  3  degree  angle  of  attack.  For  some  require¬ 
ments,  eg.,  maneuvering,  normal  forces  (such  as  lift)  may  be 
of  primary  concern  and  errors  in  drag  can  be  ignored,  but  in 
other  cases  (such  as  in  commercial  aircraft)  an  accurate  pre¬ 
diction  of  drag  is  of  paramount  concern. 

One  important  aspect  of  the  flows  around  naval  surface 
ships  and  submarines  is  the  extensive  regions  of  highly  three- 
dimensional  flow  and  especially  three-dimensional  separa¬ 
tions.  These  3-D  separations  are  gradual  developments  of 
shear  layers  which  lift  off  from  the  surface  and  wrap  into  sheets 
of  vorticity.  They  are  sensitive  to  surface  (skin)  friction  and 
thus  to  turbulence  in  the  boundary  layer.  Predictions  are  both 
important  and  difficult.  An  example  which  shows  both  the 
advances  which  have  been  made  in  the  state  of  the  art  and  the 
need  for  further  developments  is  shown  in  Fig.  1^.  The  impor¬ 
tant  qualitative  features  of  the  separations  have  been  captured, 
and  even  the  quantitative  agreement  is  impressive  over  much 
of  the  body,  but  the  pressures  on  the  important  regions  near  the 
stem  and  near  the  vortex  sheets  require  further  developments. 
The  authors  point  in  particular  to  the  need  for  further  work  in 
the  turbulence  modeling. 

Other  requirements  may  demand  accurate  prediction  of 
flow  details  such  as  the  junction  vortices  around  control  sur¬ 
faces  or  vanes  in  a  compressor.  Some  recent  work  on  tip  vortex 
flow  concluded  that  a  specialized  one-equation  turbulence 
model  worked  well  for  the  velocity  distributions,  which  are  of 
primary  importance  for  aircraft  applications.  However,  the 
pressure  in  the  vortex  core  was  severely  underpredicted  (25%) 
which  would  have  a  major  impact  on  cavitation  prediction  for 
naval  applications.  The  authors  argue  that  the  error  results 
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from  “...a  poor  representation  of  the  normal  Reynolds  stresses 
typical  of  one-equation  turbulence  models.” 

Most  of  the  examples  mentioned  so  far  have  utilized 
simple,  zero-  or  one-equation  models  for  the  turbulence  which 
we  know  require  ad  hoc  specification  of  at  least  one  charac¬ 
teristic  scale.  Two  equation  models,  the  minimum  in  complex¬ 
ity  which  do  not  require  ad  hoc  specification  of  turbulence 
scales  at  the  most  basic  level,  require  a  significant  increase  in 
computer  resources  and  care  with  numerical  stability  and 
accuracy.  Their  inherent  capabilities,  however,  should  justify 
their  use  in  practical  problems  containing  significant  flow 
complexity.  Even  here,  though,  care  must  be  taken  to  include 
sufficient  capability  demanded  by  the  problem  to  be  ad¬ 
dressed.  Fig.  2.  compares  the  development  of  turbulent  kinetic 
energy  in  homogeneous  turbulence  subject  to  suddenly  im¬ 
posed  homogeneous  shear  at  time  zero.  Standard  models, 
including  standard  two-equation  models,  are  far  from  ade¬ 
quate  in  predicting  this  flow.  A  two-equation  model  which 
includes  a  special  modification  based  on  a  perturbation  analy¬ 
sis  is  required  to  capture  the  development  of  the  flow .  It  should 
be  noted  that  the  models  are  evaluated  against  a  direct  numeri¬ 
cal  simulation  of  the  flow.  Such  comparisons  are  now  much 
more  commonplace  and  provide  needed  details  for  model 
evaluation. 

“Standard”  two-equation  models  which  do  not  incorpo¬ 
rate  strain-dependent  coefficients  cannot,  even  in  principle, 
predict  the  effects  of  streamline  curvature,  system  rotation,  or 
secondary  flows  arising  from  normal  Reynolds  stress  anisot¬ 
ropies.  In  the  past  these  deficiencies  were  overcome  by  ad  hoc 
modifications,  but  today  more  sophisticated,  but  still  two- 
equation,  models  are  being  developed  to  address  these  prob¬ 
lems.  One  difficulty,  however,  is  that  this  extensive  research 
“...has  led  to  large  numbers  of  specialized  CFD  codes  limited 
to  specific  problems”  and  thus  “...seriously  impede  the  accep¬ 
tance  of  CFD  in  the  design  community.”*^ 

Do  even  more  complicated  models  (second  order  clo¬ 
sures)  perform  better  still?  Here  at  best  we  find  confusion.  A 
recent  study'^  compared  an  advanced  two-equation  model 
with  a  standard  full  second  order  closure  (SOC)  for  a  rather 
simple  but  challenging  two-dimensional  diffuser  with  flow 
separation  at  the  wall.  The  two-equation  model  performed 
well,  but  the  more  complex  model  failed  to  predict  even  the 
occurance  of  separation!  However,  the  state  of  two-equation 
modeling  is  still  confusing  for  application  purposes;  recently 
an  engineering  manager  of  a  major  aerospace  firm  stated  that 
they  have  for  the  present  completely  “abandoned”  two-equa¬ 
tion  models  for  lack  of  confidence  in  them. 

Considering  the  hesitation  to  embrace  even  two-equation 
models,  full  SOC  is  clearly  merely  a  research  topic  at  this 
point.  The  state  of  the  art  in  this  area  has  been  summarized 
thus:***  “...existing  SOC’s  are  not  capable  of  describing  flows 
that  are  far  from  equilibrium  and  have  major  problems  with 
wall-bounded  flows.”  The  latter  point  concerning  wall- 


bounded  flows  deserves  further  discussion  because  of  its 
importance.  It  has  been  stated'^  that  the  “Achilles  heel  of 
second-order  closures  is  wall-bounded  turbulent  flows.”  The 
reasons  are  rather  arcane,  but  they  necessitate  in  many  cases 
specifying  ad  hoc  models  of  higher  order  terms  near  walls 
which  depend  on  the  distance  from  the  wall  -  a  difficulty  in 
complex  geometries.  The  author  concludes,  “Entirely  new 
non-equilibrium  models  are  needed...”  for  the  terms  causing 
the  difficulties.  It  is  further  pointed  out,  though,  that  a  SOC, 
not  a  two-equation  model,  will  be  needed  for  problems  which 
include  relaxation  or  extensive  nonlocal  effects. 


Futuire  Directions 

The  two  opposing  desires  -  CFD  codes  which  are  more 
accurate  and  general,  and  CFD  codes  which  are  faster  and 
cheaper  —  will  continue  to  define  the  directions  of  develop¬ 
ment  in  numerics,  gridding,  and  turbulence  modeling.  For 
example,  in  numerics  a  tradeoff  may  be  between  a  higher  order 
method  and  a  faster  code  requiring  less  memory.  In  gridding 
the  simplicity  of  structured  and  fixed  grids  may  compete  with 
the  generality  and  potential  overall  accuracy  of  unstructured 
and  adaptive  grids. 


Figure  2. 

Evolution  of  kinetic  energy  in  homogeneous  turbulence  sub¬ 
jected  to  uniform  shear  atf  =  0. 
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In  turbulence  modeling  the  tradeoffs  themselves  are  not 
yet  well-defined.  It  is  generally  agreed  that  higher  order  mod¬ 
els  should  permit  inclusion  of  more  physics  in  the  computa¬ 
tions.  But  the  performance  today  of  the  most  sophisticated 
models  is  mixed  at  best.  The  simpler  models  have  been  shown 
repeatedly  to  be  applicable  only  to  limited  types  of  flow  and 
not  accurate  enough  even  then  for  many  requirements. 

Developments  continue  to  occur,  however,  and  there  is 
some  hope  in  recently  proposed  models  at  the  intermediate, 
two-equation  level.  The  newer  models^  ^  are  formally  deriv¬ 
able  from  full  SOCs  and  contain  increased  generality  via 
strain-dependent  terms  and  coefficients  to  overcome  some 
deficiencies  in  earlier  models.  However,  they  are  still  strictly 
applicable  only  to  turbulence  which  is  locally  homogeneous 
and  in  equilibrium,  or  very  nearly  so. 

An  alternative,  or  perhaps  parallel,  approach  is  a  catalog¬ 
ing  of  existing  lower  order  models  and  the  applications  for 
which  they  have  been  validated.  This  is  at  best  a  temporary 
measure  unless  the  concept  of  “zonal”  modeling  is  developed, 
wherein  turbulence  models  are  selected  dynamically  depend¬ 
ing  on  the  type  of  flow  encountered.  In  the  near  term  categories 
of  models  selected  for  specific  design  problems  may  be  useful, 
but  significant  design  efforts  usually  involve  new  configura¬ 
tions  significantly  different  from  previous  ones  which  thus 
preclude  this  approach. 

In  the  intermediate  and  far  term  the  significant  advantages 
of  full  SOC  still  appear  promising.  There  are  definite  problems 
to  be  overcome,  but  formal  approximation  techniques  con¬ 
tinue  to  be  developed  and  applied  for  generating  appropriate 
models  of  higher  order  quantities.  Among  the  many  benefits 
of  the  development  of  direct  numerical  simulation  of  turbu¬ 
lence  (DNS)  over  the  past  dozen  years  or  so  is  the  ability  to 
examine  turbulence  quantities  that  have  been  impossible  to 
measure.  Such  simulations  have  led  to  important  develop¬ 
ments  in  modeling  the  higher  order  terms  in  the  SOCs.  If  the 
critical  problem  of  calculating  near  a  wall  is  solved,  the  SOC 
will  become  impossible  to  ignore,  even  though  a  major  com¬ 
puter  resource  cost  penalty  may  be  incurred. 

The  proper  representation  of  turbulence  is  indeed  a  sig¬ 
nificant  roadblock  to  the  ultimate  promises  of  CFD.  However, 
advances  on  many  fronts  -  from  cataloged,  validated  lower 
order  models  to  rationally  constructed  two-equation  models  to 
advances  in  full  second  order  closures  -  indicate  important 
new  capabilities  ahead  for  the  many  areas  of  engineering 
design  involving  the  flow  of  fluids.  The  history  of  progress  in 
this  area  is  a  prime  example  of  researchers,  developers,  and 
designers  interacting,  though  not  without  a  struggle,  to  realize 
the  promises  of  CFD. 


Biography 

Dr.  L.  Patrick  Purtell  graduated  from  the  University  of 
Maryland  with  a  Doctorate  of  Philosophy  in  Mechanical  Engi¬ 
neering.  He  is  a  Program  Officer  for  the  Office  of  Naval  Research 
and  previously  held  positions  at  the  David  Taylor  Model  Basin 
and  at  the  National  Institute  of  Standards  and  Technology. 

References 

1 .  Rouse,  H.;  Ince,  S.,  History  of  Hydraulics,  Iowa  Institute 
of  Hydraulic  Research,  1957. 

2.  Frisch,  U.  &  Orszag,  S.A.,  Physics  Today  P.  24,  January, 
1990. 

3.  Reynolds,  O.,  Phil.  Trans.  Royal  Soc.  London,  Series  A. 
Vol.  186, 1895. 

4.  Hanjalic,  K.,  Int.  J.  Heat  and  Fluid  Flow,  vol.  15,  no.  3, 
June,  1994. 

5.  Ghosal,  S.,  &  Moin,  R,  J.  comp.  Phys.,  vol.  118, 1995. 

6.  Deng,  G.B.;  Queutey,  R;  Visonneau,  M.,  Proc.  Numerical 
Ship  Hydrodynamics,  1993. 

7.  Hirschel,  E.H.;  Stock,  H.W.;  Cousteix,  J.,  Proc,  2nd  Inti. 
Symp.  Engr.  Turb.  Modeling  and  Measurements,  W.  Rodi 
&  F.  Martelli  editors,  Elsevier  Science  Publ.,  1993. 

8.  Marconi,  F.;  Sislari  M.;  Carpenter,  G.;  Chow,  R.,  AIAA- 
94-2237,  1994. 

9.  Taylor,  L.K.;  Arabshahi,  A.;  Whitfield,  D.L.,  AIAA-95- 
0313, 1995. 

10.  Dacles-  Mariani,  J.;  Zilliac, G.G.;  Chow,  J.S.;  Bradwhaw, 
R,  to  be  published,  1995. 

11.  Speziale,  C.  &  Xu,  X.-H.,  Towards  the  Development  of 
Second-order  Closure  Models  for  Non-equilibrium  Tur¬ 
bulent  Flows,  to  be  presented  at  the  Tenth  Symp.  on  Turb. 
Shear  Flows,  University  Park,  PA,  Aug.,  1995. 

12.  Melnik, R.E.;  Siclari,  M.J.;  Barber, T.;  Verhoff,  A.,  AIAA- 
94-2235, 1994. 

13.  Durbin,  P. A.,  Center  for  Turb.  Res.  Manuscript  152,  Stan¬ 
ford  Univ.,  1994. 

14.  Speziale,  C.G.,  “Modeling  of  Turbulent  Transport  Equa¬ 
tions”,  to  be  published,  1995. 

15.  Ibid. 

16.  Speziale,  C.G.,  Proc.  20th  Symp.  Naval  Hydrodynamics, 
Natl.  Acad.  Press,  to  be  published,  1995, 


Twoll995  25 


Computational  Analysis 
of  a  Missile  at  a  High 
Angle  of  Attack 
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Introduction 

Air-to-air  missiles  currently  deployed  in  the  fleet  are  not 
designed  for  high  angle  of  attack  launches  from  maneuvering 
aircraft.  The  missile  airframes  are  aerodynamically  stable  and 
“tip-off’  or  “weather-cock”  into  the  free  stream  when 
launched  at  angles  of  attack  greater  than  30  degrees.  Survival 
in  air  combat  depends  in  part  on  the  design  of  agile  missile 
airframes.  Airframe  control  must  be  maintained  even  in  very 
high  angle  of  attack  situation — aerodynamic  conditions  well 
beyond  current  design  code  capabilities. 

In  the  past,  missiles  were  designed  with  the  use  of  empiri¬ 
cal  databases,  flight  test  data,  wind  tunnel  testing,  engineering 
judgment,  and  computer  codes  based  on  panel  methods  with 
capabilities  restricted  to  moderate  angles  of  attack.  Recent 
advances  in  numerical  analysis  methods  and  advanced  com¬ 
puting  power  allow  solutions  of  three-dimensional  flowfields 
with  better  models  for  the  physics  of  the  flow.  Detailed  flow- 
field  analysis  is  provided  by  computational  fluid  dynamics 
(CFD)  methods  which  can  predict  vortex  flow  development, 
vortex  bursting,  and  flow  asymmetries  for  high  angle  of  attack 
flight  performance.  However,  these  CFD  codes  must  be  vali¬ 
dated  by  comparing  the  CFD  results  with  extensive  and  de¬ 
tailed  experimental  data  including  flowfield  information,  in 
addition  to  the  usual  surface  and  integrated  forces  and  moment 
measurements. 

At  the  higher  angles  of  attack,  which  occur  at  subsonic 
speeds,  the  flowfield  is  complex  and  dominated  by  strong 


vortices.  Viscous  effects  significandy  affect  the  growth  and 
dissipation  of  separated  and  rotational  flows,  requiring  the  use 
of  Navier-Stokes  equations  to  accurately  predict  missile  per¬ 
formance.  CFD  methods  using  the  Navier-Stokes  equations 
represent  the  most  concise  formulation  of  the  conservation 
laws  with  all  the  relevant  physics  to  predict  the  performance 
of  missile  airframes  at  high  angles  of  attack.  The  analysis 
method  also  gives  a  better  understanding  of  the  stability  and 
control  required  in  a  missile  to  safely  and  successfully  deploy 
it  from  aircraft  with  supermaneuvering  capabilities,  up  to 
70-degree  angles  of  attack. 

The  Navy  has  long  since  been  interested  in  resolving  the 
high  angle  of  attack  problem;  thus,  work  that  focuses  on 
predicting  steady  and  unsteady,  post-stall  performance  of  mar¬ 
ginally  stable  missile  airframes,  operating  at  both  high  angles 
of  attack  and  at  high  acceleration  levels,  has  begun. 

An  Analytical  Approach  to  the 
High  Angle  of  Attack  Problem 

For  the  steady  angle  of  attack  results  described  herein,  the 
OVERFLOW  code  developed  at  NASA-Ames  is  used  for  the 
numerical  simulation.  In  this  code,  the  overset  Chimera  tech¬ 
nique  is  implemented.  The  vortical  flows  are  computed  by 
using  the  Baldwin-Lomax  eddy  viscosity  model  with  the 
modifications  suggested  by  Degani  and  Schiff.* 

The  complex  surface  geometry  is  subdivided  into  simpler 
components.  The  numerical  mesh  for  the  components  such  as 


26  Naval  Research  Reviews 


Figure  1. 

Grid  Partition  for  the  Complete  Missile  Configuration. 
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fins,  tails,  and  the  missile  body  are  generated  first,  then 
brought  together  to  form  a  unit.  Partitioning  helps  reduce  the 
in-core  memory  requirements  and  allows  local  grid  refinement 
in  selected  flow  regions.  The  grid  partitioning  for  a  complete 
missile  configuration  is  shown  in  Figure  1, 

The  PEGSUS  code  was  used  to  preprocess  the  grid  data, 
providing  interpolation  stencils  at  the  grid  interfaces  and  de¬ 
fining  regions  where  holes  are  cut  by  one  grid  into  another. 
The  overset  method  provides  flexibility  in  decomposing  the 
computational  domain  and  facilitates  the  grid  generation  over 
complex  geometries,  as  well  as  small  modifications  to  the 
original  configuration.  Symmetric  solutions  for  a  half  body 
were  obtained,  using  reflection  at  the  symmetry  plane. 

The  code  uses  the  thin-layer  compressible  Navier-Stokes 
equations  to  obtain  the  numerical  solutions.  The  strong  conser¬ 
vation  law  form  of  the  governing  equations  for  a  curvilinear 
coordinate  system  (^,  T|,  Q  along  the  axial,  circumferential,  and 
normal  direction,  respectively,  is  as  follows: 


Figure  2. 

Missile  Configuration  With  Full  Wings. 
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3,Q  +  a|F  +  9^G  +  9^H  =  Re'8^S  (1) 

The  numerical  integration  is  performed  using  a  partially 
flux-split  numerical  scheme.  A  two-factor  scheme  is  selected 
which  uses  Steger- Warming  flux  vector  splitting  for  the  stream 
direction  and  central  differences  for  the  other  two  directions. 
For  the  central  differencing  directions,  implicit  and  explicit 
smoothing  are  used.  The  resulting  two-factored  algorithm  is 

[{I  +  h?^  (Ay  +  h^C"  -  hRer\r'M"J  - 
x[I  +  h^^  (A-  r  +  -  Dii^]  Aq" 

=  -  At  {  8^  [n"  -  F"-]  +  8^  -  FI 

+  8^(G"-G^)  +  8^^/F-//„j 

+  (2) 

In  Eq.  (2),  De  are  the  explicit  dissipation  terms  that  are 
used  along  the  directions  where  central  differencing  is  em¬ 
ployed,  and  Di  are  the  implicit  dissipation  terms  that  are  added 
for  numerical  stability.  The  dissipation  terms  used  are  a  com¬ 
bination  of  second-  and  fourth-order  terms.^ 

Steady  CFD  Results  for  an 
Advanced  Configuration 

Using  the  analytical  approach  described  above,  the  mis¬ 
sile  flowfield  for  an  advanced  missile  configuration  was  com¬ 
puted  for  45-degrees  roll,  which  is  the  preferred  orientation 
during  flight  requiring  maneuvers  to  the  target.^  The  missile 
configuration  considered  is  shown  in  Figure  2.  Solutions  were 
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Figuro  3. 

Missile  Forebody  and  Fin  Gap  Grids. 


obtained  for  a  free-stream  Mach  number  of  0.3,  angles  of 
attack  up  to  60  degrees,  and  a  Reynolds  number  of  1  million, 
which  is  typical  of  a  high  angle  of  attack  launching  scenario. 
Flowfield  solutions  were  computed  with  the  fins  off  and  on. 
The  gap  between  the  body  and  fin  was  also  modeled,  as  shown 
in  Figure  3.  The  Reynolds  number,  based  on  the  missile 
diameter,  was  high  enough  to  assume  fully  turbulent  flow 
which  was  used  in  the  analysis. 

Shaded  contours  of  Mach  numbers  are  shown  in  Figure  4 
at  a  crossflow  plane  aft  the  forward  canards.  The  flow  is 
coming  from  the  lower  right  to  upper  left  in  the  figure,  with 
the  missile  body  at  a  45-degree  angle  of  attack  and  both 
canards  deflected  an  additional  10  degrees.  The  darker  blue 
regions  above  the  missile  body  are  the  result  of  flow  separation 
over  the  upper  and  lower  canards.  However,  where  the  fins  do 
not  influence  the  flow  on  the  windward  side  of  the  body,  an 
accelerated  flow  region  is  shown  in  red.  The  characteristics  of 
shear  flow  are  shown  in  the  small  area  between  the  two  flow 
regions. 

For  the  same  flow  conditions,  the  vortex  system  with 
helicity  (measure  of  the  vortex  strength)  densities  at  several 
crossflow  planes  is  shown  in  Figure  5.  Flow  over  the  fins  was 
shown  to  be  highly  separated  with  no  coherent  leading  edge 
vortex.  The  forebody  vortex,  as  well  as  the  vortices  generated 
in  the  gap  between  the  fins  and  the  body,  is  complicated, 
merging  into  a  single  vortical  flowfield  region  over  the  upper 
fin.  Downstream  of  the  canards,  vortices  are  generated  by  the 
missile  body  which  did  not  lift  away  from  the  missile  surface, 
but,  rather,  merged  with  the  tail  vortical  flowfield. 

Normal  force  and  pitching  moment  coefficients  for  ex¬ 
perimental  data,  computations,  and  Missile  DATCOM  (Data 
Compendium  Air  Force  generated  aerodynamic  database)  are 
shown  in  Figure  6  for  the  missile  in  the  “x”  orientation  with 
canard  fins  deflected  10  degrees.  For  this  case,  the  1993 
Missile  DATCOM  database  gives  good  results;  however,  no 


Figure  4. 

Mach  Number  Contours  in  Crossflow  Plane  Aft  of  Canards. 
M.  =  O.S.a  =  45°.  (()  =  45°,  Scan  =  1 0°,  Reo  =  .95  x  1 0®. 


knowledge  is  gained  about  the  details  of  the  flowfield  to  assist 
in  making  rapid  design  improvements. 

The  computed  flowfield  over  the  missile  indicates  that  the 
complex  flow  pattern  generated  by  the  deflected  fins  plays  an 
important  role  in  the  development  of  the  vortical  flowfield 
over  a  missile  at  high  angles  of  attack.  In  the  future,  complete 
body  modeling  may  be  required  because  of  the  asymmetric 
vortices  that  are  known  to  be  generated  by  long  slender  bodies 
at  higher  angles  of  attack  to  get  better  agreement  with  experi¬ 
mental  data. 

For  the  missile  in  the  “+”  orientation  with  zero  canard  fin 
deflection,  shown  in  Figure  7,  computations  executed  with  the 
General  Aerodynamic  Simulation  Program®  (GASP)  give  ex¬ 
cellent  agreement  with  experiments,  while  Missile  DATCOM 
grossly  underpredicts  both  the  normal  force  and  pitching 
moment  coefficients.''  The  details  of  GASP  calculations  are 
beyond  the  scope  of  this  discussion,  therefore  are  not  included 
herein. 

Unsteady  CFD  Results  for 
Pitching  Missile  Body 

ANavier-Stokes  simulation  of  laminar,  low  subsonic  flow 
over  a  pitching  missile  body  was  executed  for  incidence  angles 
up  to  25  degrees.  The  CFL3D®  code  was  used  to  obtain  the 
numerical  solutions.  (Details  on  this  code  can  be  found  in 
Reference  5).  The  flow  was  assumed  to  be  symmetric  with 
respect  to  the  pitching  plane  throughout  the  calculations. 
Computed  results  of  helicity  at  several  crossflow  planes  along 
the  body  are  given  for  a  missile  body  pitching  at  a  high  rate  of 
500  deg/s.® 

Figure  8  gives  the  computed  helicity  for  a  missile  body 
alone  at  the  peak  angle  in  the  oscillation  cycle  shown  in  the 
shape  of  a  horseshoe  around  the  body.  A  prominent  vortical 
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Figure  5. 

Vortical  Flow  Field  over  the  Missile  Shown  by  Helicity 
Ai.  =  0.3, a  =  45°, «])  =  45°,  dean  =  10°,  Reo  =  .95  x  10°. 


Fin  Vorttces 


flow  structure  is  shown  at  a  25-degree  angle  of  attack  with 
vortices  beginning  to  spread  away  from  the  body  because  the 
pitch  rate  has  decreased  to  0  deg/s  at  the  top  of  the  oscillation 
cycle. 

The  computations  clearly  show  the  vortex  lag  effects 
during  the  oscillation  cycle,  as  shown  in  Figures  9  and  10, 
where  the  vortex  strengths  are  larger  on  the  downstroke 
when  compared  with  those  at  the  same  angle  of  attack  on 
the  upstroke.  However,  the  missile  body  produces  more  lift 
on  the  upstroke  at  15  degrees  when  the  vortices  are  tightly 
concentrated  around  the  leeward  side  of  the  body,  as  op¬ 
posed  to  the  large  circular  and  spread  vortical  flow  patterns 
for  15  degrees  on  the  downstroke,  as  shown  in  Figure  9. 
Only  the  start  of  vortical  flow  development  at  the  last 
crossflow  location  is  shown  in  Figure  10  at  a  5-degree  angle 
of  attack  on  the  upstroke,  but  vortical  flow  still  exists  at  all 
crossflow  locations  for  5-degree  angle  of  attack  on  the 
downstroke. 

An  extension  of  the  unsteady-body-alone  work  is  in  pro¬ 
gress  for  the  advanced  missile  configuration  with  undeflected 
canards  and  tail  fins  in  the  same  angle  of  attack  range  and  at  a 
more  realistic  pitching  rate  of  250  deg/s. 

A  Look  at  the  Future 

Two-dimensional  viscous  analyses  with  Reynolds-aver¬ 
aged  Navier-Stokes  equations  are  easily  computed  and  give 
reliable  results  for  most  aerodynamic  design  purposes.  Over 
the  years,  it  has  become  evident  that  numerical  solutions, 
including  a  transition  model,  most  closely  match  the  extensive 


experimental  results  for  a  fin  at  high  angles  of  attack  in  the 
transitional  Reynolds  number  flight  regime. 

In  recent  years,  three-dimensional  viscous  analysis  has 
become  a  more  common  occurrence,  enabled  by  more  ad¬ 
vanced  computational  vector  processing.  In  fact,  numerical 
investigations  have  been  performed  to  compute  the  flow 
over  delta  and  double-delta  wings  at  high  angles  of  attack, 
and  successful  efforts  have  been  conducted  to  obtain  major 
flow  features  with  reasonable  accuracy,  including  the  very 
difficult  phenomenon  of  vortex  breakdown.  In  addition, 
preliminary  computations  have  been  successfully  com¬ 
pleted  for  an  advanced  missile  configuration  at  several 
steady  angles  of  attack  using  the  OVERFLOW  code  with 
various  turbulence  models. 


Figure  6. 

Comparison  of  Measured  and  Computed  Normal  Forces  and 
Pitching  Moment  Coefficients  ("x"  orientation). 
U.  =  0,3, <))  =  45°,  dean  =  10°,  Reo  =  .95  x  10°, 
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Figure  7. 

Comparison  of  Measured  and  Computed  Normal  Forces  and 
Pitching  Moment  Coefficients  orientation). 
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As  CFD  capabilities  advance,  so  does  the  confidence  in 
this  tool.  Where  in  the  past,  wind  tunnel  tests,  flight  tests,  and 
simulations  were  the  workhorses  in  the  development  of  mis¬ 
siles,  today  it  is  becoming  more  cost  effective  and  advanta¬ 
geous  to  rely  more  on  CFD  work  and  less  on  the  other  methods, 
although  some  of  each  is  still  necessary. 

As  a  direct  result  of  CFD  advances,  design  capability  is 
quickly  improving.  In  the  near  future,  it  will  allow  designers 
to  rapidly  predict  performance  and  conduct  static  or  dynamic 
stability  analysis  in  the  various  design  stages  in  an  efficient 
and  reasonably  accurate  manner.  Understanding  the  mecha¬ 
nisms  that  generate  vortex  breakdown  and  the  ability  to  predict 
its  occurrence,  for  example,  are  crucial  to  efforts  for  improving 
missile  designs  and  missile  performance  in  flight  at  high 
angles  of  attack. 
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Computed  Helicity  at  Selected  Missile  Body  Crossflow  Planes, 
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Introduction 

Combustion  is  a  complex  phenomenon  involving  fluid 
dynamic  mixing  between  the  fuel  and  oxidizer,  chemical  reac¬ 
tions  with  the  associated  kinetics  and  thermodynamics,  and 
heat  transfer  to  the  system  and  surroundings.  It  is  further  compli¬ 
cated  by  the  complex  geometries  of  the  combustion  chambers  in 
use  today.  Current  and  future  applications  require  increased 
energy  release  with  reduced  chamber  volume,  increased  equilib¬ 
rium  temperatures,  multiphase  reacting  flows  with  radiative  heat 
transfer  and  sometimes  even  electric  and  magnetic  fields.  Basic 
understanding  of  the  physical  and  chemical  processes  involved 
in  combustion  becomes  a  necessity  for  the  design  of  efficient  and 
reliable  combustion  systems. 

Traditional  development  of  combustion  systems  based 
largely  on  empirical  data,  and  relatively  simple  design  models 
will  not  be  sufficient  if  energy  utilization  is  to  be  maximized. 
Detailed  full-scale  experiments  have  become  exceedingly  expen¬ 
sive  and  often  times  very  difficult.  This  is  particularly  so  when  a 
parametric  evaluation  is  desired.  On  the  other  hand,  "numerical 
experiments"  open  an  avenue  for  isolating  and  understanding  the 
influence  of  different  parameters  in  complex  combustion  situ¬ 
ations.  It  also  provides  a  potential  design  tool  for  future  combus¬ 
tion  and  propulsion  systems.  However,  this  is  not  simple  or 
straight  forward  as  will  be  discussed  below. 


Rapid  advances  in  super  computing  and  massively  paral¬ 
lel  processing  has  brought  forth  computational  capabilities 
hitherto  unheard  of.  Successful  predictions  of  the  detailed  flow 
fields  in  complex  aerodynamic  systems  have  been  demon¬ 
strated,  and  employed  in  design.  Extension  of  Computational 
Fluid  Dynamics  (CFD)  to  combustion,  oftentimes,  has  been 
considered  as  adding  a  few  more  terms  in  the  governing  fluid 
flow  equations.  However,  when  one  realizes  that  combustion 
involves  utilization  of  fuels  and  other  species  of  vastly  differ¬ 
ent  molecular  properties,  chemical  reactions  with  vasdy  dif¬ 
ferent  time  scales,  wide  range  of  stoichiometry  and  load 
conditions,  and  slurries,  gels,  solids  and  gas-particle  systems 
with  precise  thermal  and  exhaust  emission  management  (for 
pollution  and  signature  control),  the  complexity  is  obvious. 

Because  of  the  complexity  of  the  problem,  it  is  not  feasible 
to  directly  simulate  the  turbulent  reacting  flow  in  a  typical 
Navy  propulsion  system.  Simplifying  assumptions,  some¬ 
times  drastic,  need  to  be  made.  The  simplifications  will  vary 
with  the  problem  and  in  an  industrial  context  may  also  depend 
on  economic  and  temporal  constraints.  The  resources  allocated 
for  computations  and  turn-around  time  needed  is  quite  differ¬ 
ent  in  a  design  process  than  in  a  point  design  analysis  or 
failure-analysis  scenario.  The  level  of  detail  included  in  the 
numerical  models  and  the  simplifying  assumptions  made  will 
also  be  different  depending  on  the  problem  to  be  solved  and 
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the  level  of  detail  expected  in  the  results.  Therefore,  in  this 
paper,  instead  of  focussing  on  a  specific  scenario  of  perhaps 
greater  practical  interest,  we  address  the  general  problem  of 
computational  combustion.  By  doing  so,  we  hope  the  reader 
will  have  an  appreciation  of  the  complexity  of  the  combustion 
phenomena  and  the  need  to  carefully  choose  the  simplifying 
assumptions  needed  in  particular  applications. 

The  Office  of  Naval  Research  (ONR)  has  been  sponsoring 
studies  on  computational  combustion  involving  gaseous,  liq¬ 
uid  and  solid  fuels,  and  presently  a  new  class  of  strained 
high-energy  hydrocarbon  fuels.  Substantial  accomplishments 
have  been  made,  and  the  scientific  understanding  of  combust¬ 
ing  flows  has  been  greatly  enhanced  by  these  studies.  An 
attempt  is  made,  in  this  paper,  to  address  the  various  issues 
involved  and  the  progress  made  in  computationally  simulating 
combustion  problems  of  relevance  to  the  Navy,  with  an  em¬ 
phasis  on  research  sponsored  by  ONR, 

For  most  Navy  propulsion  applications,  combustion  takes 
place  in  a  turbulent  flow.  Hence,  the  approaches  adopted  to 
compute  a  turbulent  flow,  discussed  elsewhere  in  this  issue, 
are  of  interest  to  turbulent  combustion.  However,  even  for  low 
Mach  number  applications,  incompressible  flow  assumptions 
are  not  appropriate  for  combustion  because  of  the  large  densi  ty 
changes  and  pressure  waves  generated  by  the  energy  released. 
In  this  report,  we  first  discuss  the  various  physical  and  chemi¬ 
cal  processes  that  must  be  considered  in  a  general  combustion 
problem.  This  brings  out  the  complexity  of  the  combustion 
phenomenon  and  the  difficulties  in  trying  to  solve  directly  the 
governing  equations.  A  laminar  flame  problem  is  used  as  an 
illustrative  example  since  it  one  of  the  problems  that  can  be 
solved  directly  and  is  often  used  to  evaluate  and  calibrate 
simplified  models  for  use  in  more  practical  Navy  problems. 
This  is  followed  by  brief  discussions  of  various  approaches 
adopted  in  turbulent  combustion  in  work  sponsored  by  ONR 
and  the  progress  made  recently.  Irrespective  of  the  approach 
used  for  a  particular  problem,  whether  it  be  modeling  or 
simulation,  issues  such  as  chemistry  solvers  and  models,  mul¬ 
tiphase  flow  effects,  soot,  and  heat  transfer  including  radiation 
need  to  be  considered  in  practical  combustion  problems. 
Therefore,  the  progress  made  in  these  areas  are  also  high¬ 
lighted  with  specific  examples.  Future  prospects  and  research 
needs  in  computational  combustion  are  discussed  throughout 
and  highlighted  in  the  final  section. 

Computational  Combustion 

Computational  combustion  is  sometimes  thought  of  CFD 
with  an  additional  source  term  in  the  energy  equation  to 
account  for  the  heat  released  in  chemical  reactions.  Although 
the  overall  exothermicity  is  what  makes  combustion  of  interest 
to  Navy  propulsion,  there  are  various  other  physical  and 
chemical  processes  involved  in  combustion  that  must  be  un¬ 
derstood  to  optimize  and  control  the  energy  release  in  order  to 


safely,  cleanly  and  efficiently  operate  the  propulsion  systems. 
Some  of  the  physical  and  chemical  processes  that  must  be 
considered  in  a  typical  combustion  system  is  illustrated  sche¬ 
matically  in  Figure  1. 

The  energetic  decomposition  of  a  fuel  or  the  reaction  of  a 
fuel  with  an  oxidizer  produces  a  large  number  of  intermediate 
species  which  need  to  be  identified  and  tracked.  The  properties 
of  these  individual  species  such  as  their  molecular  weights  and 
chemical  properties  are  quite  distinct.  Therefore,  as  many 
convective  transport  equations  as  there  are  species  are  needed. 
In  addition,  these  species  tend  to  diffuse  into  each  other, 
conduct  heat  to  varying  degrees  and  some  of  them  may  also 
radiate.  Therefore,  diffusive  transport  processes  such  as  multi¬ 
species  diffusion,  heat  conduction,  viscous  and  thermal  diffu¬ 
sion  need  to  be  considered.  Furthermore,  all  these  processes 
are  coupled  together  and  may  interact  with  each  other  non 
linearly.  In  general,  combustion  takes  place  in  confined  ge¬ 
ometries  for  propulsion  applications,  and  therefore,  heat  trans¬ 
fer  to  the  walls,  boundary  layers  and  surface  reactions  must 
also  be  considered.  Quite  often,  the  fuel  is  stored  in  the  form 
of  a  liquid.  Therefore,  the  atomization,  droplet-gas  interactions 
and  vaporization  of  fuels  need  to  be  accounted  for.  Depending 
on  the  specific  problem,  external  forces  on  the  system  such  as 
buoyancy  forces  and  magnetic  and  electrical  fields  may  also 
need  to  be  included.  Because  of  the  high  temperatures  and 
pressures  encountered  in  some  applications,  various  equations 
of  state  may  also  need  to  be  considered. 

The  complexity  of  the  problem  of  computational  combus¬ 
tion  is  what  also  provides  one  of  its  major  strengths.  Since 
different  interacting  processes  are  accounted  for,  we  can  gain 
better  understanding  of  many  combustion  problems  by  turning 
processes  on  and  off  or  by  assigning  desired  values  for  them. 
For  example,  the  behavior  of  flames  or  fires  in  space  can  be 
simulated  by  assigning  low  or  a  zero  value  for  the  gravitational 
acceleration  term  in  the  governing  equations.  That  is,  a  com¬ 
putational  combustion  model  can  be  used  as  a  tool  to  perform 
numerical  experiments  which  are  difficult,  expensive  or  even 


Figure  1. 

A  schematic  illustration  of  some  of  the  physical  and  chemical 
processes  in  a  typical  combustion  system. 
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impossible  to  perform  in  the  laboratory.  These  numerical 
experiments  can  provide  valuable  information  for  the  design 
of  new  and  innovative  combustion  devices  and  to  improve 
existing  designs. 

Basic  Numerical  Requirements 

It  is  not  just  the  presence  of  additional  terms  in  the 
governing  equations  that  make  the  computation  of  combustion 
problems  more  complex  than  CFD,  but  the  numerical  require¬ 
ments  imposed  by  the  underlying  processes  represented  by 
these  additional  terms.  The  basic  numerical  requirements  dis¬ 
cussed  here  are  applicable  to  both  turbulent  and  laminar  flame 
situations  but  the  focus  in  the  discussion  below  is  more  on  a 
freely-propagating  laminar  flame.  A  laminar  flame  is  of  direct 
interest  to  the  Navy  because  some  of  the  most  difficult  fires  to 
put  out  on  ships  are  those  from  small  fuel  spills.  They  also  have 
other  dual-use  applications  such  as  burners.  In  addition  to 
practical  combustion  problems,  a  more  important  use  of  the 
laminar  flame  is  as  a  unit  problem  to  test  and  evaluate  sub¬ 
models  such  as  for  chemistry  that  is  needed  in  other  practical 
applications  of  interest  to  the  Navy.  Adetailed  model  of  a  flame 
must  contain  accurate  representations  of  at  least  the  convec¬ 
tive,  diffusive,  and  chemical  processes. 

Diffusive  transport  processes  play  a  very  important  role 
in  flames  and  must  be  properly  resolved.  Numerically  this 
means  that  any  numerical  diffusion  in  the  calculation  must  be 
considerably  less  than  any  important  physical  diffusion  effect. 
A  Lagrangian  method  is  ideally  suited  for  the  convective 
transport  because  eliminating  the  advection  term  in  effect 
means  eliminating  numerical  diffusion  from  the  calculation. 
However,  extending  Lagrangian  methods  to  multidimensions 
is  both  very  difficult  and  expensive.  Therefore  Eulcrian  meth¬ 
ods  have  to  be  considered  for  multidimensional  flames. 
Eulerian  methods  have  been  successfully  used  to  simulate 
flames,  but  these  methods  need  finer  grids  than  their  La¬ 
grangian  counterparts  in  order  to  minimize  numerical  diffu¬ 
sion.  Minimizing  numerical  diffusion  is  also  important  for 
maintaining  the  sharp  gradients  present  in  flames.  Therefore, 
an  important  factor  in  the  simulation  of  flame  propagation  is 
the  choice  of  the  spatial  grid  size  or  resolution. 

Spatial  resolution  is  achieved  by  having  a  large  number 
of  grid  points  in  regions  where  the  gradients  are  steep.  In  a 
detailed  flame  model,  steep  gradients  in  different  species 
usually  occur  at  different  spatial  locations  necessitating  the  use 
of  fine  resolution  all  across  the  flame.  The  simplest  gridding 
approach  then  would  be  to  estimate  the  smallest  grid  size 
required  to  resolve  the  feature  of  the  flame  under  investigation 
and  use  this  grid  size  over  the  whole  computational  domain. 
The  cost  of  such  an  approach  is  usually  prohibitive  and  there¬ 
fore  a  non-uniform  grid  is  used  where  only  a  part  of  the 
computational  domain  is  finely  gridded  (with  small  grid  size). 
In  a  time-dependent  problem,  the  flame  may  change  its  spatial 
location  and  then  an  adaptive  gridding  approach  in  which  the 


finely  gridded  region  moves  with  some  characteristic  feature 
of  the  flame  is  chosen.  This  is  feasible  in  one-  or  two-dimen¬ 
sional  flows  but  extremely  complex  in  three-dimensional 
flows  in  complex  geometries. 

Since  flames  are  slow  (when  compared  to  the  speed  of 
sound),  one  would  like  to  use  large  time-steps  to  study  the 
dynamics.  However,  if  an  explicit  method  is  used  to  advance 
the  solution  in  time,  very  small  time  steps  which  are  restricted 
by  the  speed  of  sound  must  be  taken  to  keep  the  calculations 
numerically  stable.  Explicit  methods  are  those  in  which  the 
solution  at  the  current  time-step  depends  only  the  solutions  at 
the  previous  steps.  Explicit  methods  are  easier  to  code  and  the 
cost  per  step  is  lower  than  implicit  methods.  Implicit  methods 
are  those  in  which  the  solution  at  the  current  time-step  depends 
also  on  a  previous  estimate  of  the  solution  at  the  current 
time-step.  Thus  it  requires  an  iterative  procedure  or  matrix 
inversion  and  is  more  expensive  per  time-step.  However,  for 
slowly  evolving  problems,  implicit  methods  are  more  efficient 
because  they  require  fewer  but  larger  time-steps  to  obtain  the 
solution  for  a  given  time  interval. 

Simulating  in  detail  the  combustion  of  complex  hydrocar¬ 
bon  fuels  could  easily  involve  hundreds  of  species  and  thou¬ 
sands  of  elementary  reactions  among  them.  Solving  a  set  of 
detailed  chemical  reaction  rate  equations  is  usually  the  most 
expensive  part  of  detailed  flame  models.  In  general,  the  nu¬ 
merical  accuracy  required  of  chemistry  algorithms  is  lower  in 
a  fluid  flow  situation  than  in  a  purely  chemical  kinetics  prob¬ 
lem.  So  the  emphasis  on  chemistry  solvers  is  more  on  speed 
than  on  higher  orders  of  accuracy.  Some  of  the  recent  devel¬ 
opments  in  faster  chemistry  solvers  on  parallel  computers  are 
discussed  later. 

The  last  but  most  important  part  of  flame  models  is  the 
coupling  between  the  various  processes.  One  approach  is  to 
solve  all  the  equations  simultaneously.  This  requires  some 
form  of  iterative  method  because  the  equations  are  nonlinear. 
In  general,  this  is  a  very  expensive  approach  and  difficult  to 
efficiently  parallelize  but  it  has  been  adopted  successfully  in 
many  specific  situations  such  as  the  determination  of  burning 
velocity  and  a  steady  flame-structure.  Another  approach  is  to 
use  an  operator  or  time-step  split  procedure  in  which  the 
individual  processes  are  solved  using  the  most  efficient  meth¬ 
ods  and  then  the  solutions  are  appropriately  coupled.  Here 
particular  care  must  be  taken  in  the  coupling  procedure.  This 
approach  has  been  adopted  successfully  in  a  wide  variety  of 
reactive-flow  calculations'  and  some  examples  are  discussed 
later  in  this  paper.  A  recent  paper  compares  an  operator-split 
scheme  and  unsplit  scheme  for  the  simulation  of  one-dimen¬ 
sional  flames  and  comes  to  the  conclusion  that  the  split  ap¬ 
proach  is  slightly  superior  to  the  unsplit  approach  in  problems 
involving  detailed  chemistry^.  Various  approaches  in  between 
the  process-splitting  and  the  fully  coupled  approaches  have 
also  been  taken  and  shown  to  give  satisfactory  results  in  many 
specific  situations. 
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Figure  2. 

Comparison  of  the  instantaneous  reaction-rate  distributions 
from  a  reactive  Jet  simulation  \Arith  distributions  of  an  interme¬ 
diate  chemical  species  (OH)  used  to  characterize  the  experi¬ 
mental  combustion  regions.  (Grinstein,  Ref.  5) 
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Turbulent  Flows 

Even  with  the  rapid  advances  in  computers  and  computa¬ 
tional  methodologies,  it  is  not  currently  possible  to  simulate 
combustion  in  a  complex  propulsion  device  by  directly  solving 
the  basic  equations  discussed  above,  because  the  flow  is  usu¬ 
ally  turbulent.  The  difficulties  associated  with  and  the  issues 
that  need  to  be  addressed  in  simulating  turbulent  flows  are 
discussed  elsewhere  in  this  issue.  They  are  relevant  to  turbu¬ 
lent  combustion  and  will  not  be  elaborated  on.  Here,  attention 
is  focussed  on  additional  issues  that  are  specific  to  combustion. 

Numerical  simulations  of  turbulent  combustion  must  be 
capable  of  resolving  a  broad  range  of  spatial  and  temporal 
scales  including  multi-species  transport.  Direct  numerical 
simulations  (DNS)  of  combustion  which  involves  directly 
solving  the  basic  equations  can  be  performed  only  for  a  small 
class  of  problems,  typically  involving  low  Reynolds  (Re) 
numbers.  For  other  problems  of  interest,  simplifying  assump¬ 
tions  need  to  be  made.  One  of  the  simplifying  assumptions  that 
is  often  made  in  turbulent  combustion  is  the  assumption  of  a 
single-step  global  reaction  between  the  fuel  and  the  oxidizer. 
This  not  only  simplifies  the  chemical  kinetics  but  also  reduces 
the  number  of  species  involved  to  a  fuel,  oxidizer,  product(s) 
and  usually  a  diluent.  There  are  a  number  of  papers  dealing 
with  a  single-step  reaction  in  conjunction  with  the  DNS  of  the 
underlying  fluid  flow  in  very  simple  systems^  Research"^ 
sponsored  by  ONR  also  considers  the  temperature  and  species 
dependent  diffusive  transport  processes  and  realistic  inflow 


and  outflow  boundary  conditions.  This  paper  is  part  of  a  larger 
effort^  at  the  Naval  Research  Laboratory  (NRL)  to  use  numeri¬ 
cal  simulations  to  study  the  dynamics  and  topology  of  non-ax- 
isymmetric  jets  because  of  their  enhanced  characteristic 
entrainment  properties  relative  to  those  of  round  jets.  It  shows 
that  direct  simulations  are  invaluable  for  obtaining  a  compre¬ 
hensive  understanding  of  phenomena  such  as  axis-switching 
and  its  role  in  enhancing  the  mixing  of  non-circular  jets  over 
circular  jets. 

However,  for  higher  Re  number  unsteady  reactive  flows 
in  larger  systems,  it  also  brings  out  the  need  to  adopt  a 
large-eddy  simulation  (LES)  approach.  In  the  LES  approach, 
the  flow  is  decomposed  into  two  parts.  The  first  part  contains 
the  large  scales  and  depends  on  the  boundary  conditions  of  the 
flow.  The  second  part  contains  the  small  scales  and  is  modeled 
using  subgrid  scale  turbulence  models.  Adequate  subgrid  scale 
models  for  compressible,  multi-species  reactive  flows  do  not 
exist  now.  Therefore,  models  from  CFD^  or  those  implicit  in 
non-linear  algorithms^  such  as  the  monotonically  integrated 
large-eddy  simulation  (MILES)  are  usually  used.  In  either 
case,  comparisons  and  calibration  with  experiments  need  to  be 
performed.  For  example,  in  Fig.~2,  the  instantaneous  chemical 
reaction-rate  distributions  from  a  reactive  jet  simulation  is 
compared  with  distributions  of  an  intermediate  chemical  spe¬ 
cies  (OH)  used  to  characterize  the  experimental  combustion 
regions.  In  this  particular  example,  the  axisymmetry  of  a 
circular  jet  is  broken  up  by  introducing  vortex  generators  or 
tabs  along  the  edge  of  the  nozzle.  This  provides  an  alternate  to 
using  non-circular  nozzles  to  enhance  the  mixing  of  a  fuel  jet 
with  the  surrounding  air.  With  the  confidence  gained  by  com¬ 
parisons  to  experimental  data,  numerical  simulations  can  be 
used  to  gain  insight  into  the  details  of  mixing  process  using 
information  such  as  the  time  evolving  vorticity  field  which 
cannot  be  obtained  in  experiments. 

Classical  approaches  to  modeling  turbulence  based  on  a 
statistical  treatment  of  fluctuations  about  a  stationary  or  slowly 
varying  state  have  also  been  adapted  to  combustion^.  As  in 
modeling  non-reactive  turbulent  flows,  the  primary  variables 
in  the  equations  are  divided  into  two  parts:  a  mean,  time-aver¬ 
aged  part  of  the  flow  and  a  fluctuating  part  representing 
deviations  from  the  mean.  Models  relating  the  fluctuations  to 
the  known  quantities  of  the  flow  must  be  postulated  to  close 
the  system  of  equations.  The  closures  from  the  Reynolds 
averaged  Navier  Stokes  (RANS)  approaches  in  fluid  dynamics 
are  not  sufficient  because  of  the  additional  processes  that  must 
be  considered  in  turbulent  combustion.  The  large  changes  in 
density  that  is  typical  of  combustion  problems  also  suggests 
the  use  of  averaging  procedures  such  as  the  Favre  or  mass-av¬ 
eraged  approach.  As  part  of  the  ONR  propulsion  program, 
stochastic  modeling  of  turbulent  flows  is  also  sponsored  in 
addition  to  DNS  and  LES^.  For  example,  the  maximum  rate  of 
mean  reactant  conversion  in  a  homogeneous  turbulent  reacting 
flow  problem,  obtained  from  a  model  is  compared  to  DNS  data 
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in  Fig.-3.  This  also  brings  out  another  use  of  DNS,  as  a  tool 
to  perform  idealized  simulations  (unit  problems)  to  evaluate 
and  calibrate  models  for  use  in  more  practical  computations. 
More  recently^®,  in  this  research  effort  both  probability  density 
function  (pdO  and  moment  methods  are  considered  with  a 
focus  on  1)  mathematical  modeling  of  the  conditional  statistics 
of  scalars  in  turbulent  mixing  and  reaction  and  2)  development 
of  algebraic  closures  for  the  turbulent  flux  of  scalars  trans¬ 
ported  in  reacting  flows. 

Even  approaches  based  on  solving  the  time-averaged 
multidimensional  equations  with  phenomenological  models 
for  closure  are  often  too  expensive  (time  consuming)  and  not 
accurate  enough  for  routine  use  in  the  design  of  combustors^  ^ 
The  more  common  approach  adopted  in  the  industry  is  based 
on  semi-empirical  correlations^^  Significant  advances  still 
need  to  be  made  in  simulating  turbulent  combustion  for  it  to 
be  used  on  a  routine  basis  in  an  industrial  environment  driven 
more  by  cost  than  by  technological  advancements.  With  the 
current  emphasis  on  timed  fuel  injection  for  increased  fuel 
efficiency  and  reduced  emissions,  unsteady  aspects  of  the  flow 
will  have  to  be  resolved.  This  necessitates  the  development  of 
a  new  class  of  turbulence  models  that  are  applicable  for 
unsteady  flows. 

The  approach  adopted  to  compute  the  flow  in  a  combustor 
or  idealized  representation  of  it,  will  depend  on  the  specific 
problem  under  consideration.  Time  and  cost  constraints  will 
always  necessitate  the  use  of  models  and  approximations  in 
practical  combustion  scenarios.  In  some  cases,  the  focus  may 
be  on  combustion  instability  and  the  role  of  large-scale  vortex 
dynamics.  In  this  case,  resolving  the  inviscid  vortex-dynamics 
and  acoustics  by  solving  the  unsteady  Euler  equations  may  all 
that  be  required.  Irrespective  of  the  approach  adopted,  in 
combustion,  additional  processes  such  as  chemical  reactions 
and  heat  release,  soot  formation  and  destruction  and  heat 
transfer  must  be  represented  in  terms  of  sub-models.  Issues 
related  to  these  processes  and  recent  advances  made  in  the 
research  sponsored  by  ONR  are  discussed  next. 

Reduced  Chemistry  Models 

The  combustion  of  fuels  involve  the  decomposition  of  the 
fuel,  formation  of  intermediate  species  and  reactions  among 
them  and  the  oxidizer  to  release  chemical  energy.  Tens  of 
species  and  hundreds  of  reactions  among  them  are  needed  to 
represent  the  combustion  of  typical  hydrocarbon  fuels  used  in 
the  Navy.  In  many  situations,  useful  information  can  be  ob¬ 
tained  by  considering  the  species  to  be  in  chemical  equilibrium 
and  ignoring  the  reaction  kinetics.  There  are  many  other  situ¬ 
ations  where  the  flow  times  and  dynamics  are  such  that, 
chemical  equilibrium  does  not  exist.  However,  solving  a  de¬ 
tailed  set  of  elementary  reactions  in  conjunction  with  a  com¬ 
plex  flow  is  prohibitively  expensive.  Therefore,  there  is  a  need 
to  develop  reduced  chemistry  models. 


Figure  3. 

Comparison  of  Model  data  (-)  to  DNS  (  q  )  for  the  temporal 
variation  of  normalized  mean  fuel  concentration  in  a  homoge¬ 
neous  turbulent  reacting  flow.  (Givi  et  ai,  Ref.  9) 
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For  example,  with  the  current  emphasis  on  environmen- 
tally-sound  ships,  there  is  a  need  for  a  predictive  tool  for  the 
emission  of  NOx  from  Diesel  Engines  and  Gas  turbines.  The 
work  of  Williams’^  uses  a  four-step  model  for  fuels  of  interest 
in  Diesel  combustion.  In  this  approach,  the  first  step  is  an 
overall  fuel-consumption  step  that  produces  CO,  and  CO2. 
The  second  step  is  the  water-gas  shift  or  CO  consumption 
step.The  third  and  fourth  steps  are  associated  with  the  radical 
recombination  and  oxygen  consumption.  The  rates  of  these 
four  reduced  steps  are  derived  from  the  rates  of  the  elementary 
reactions. 

In  addition,a  one-step  description  of  the  NOx  production 
is  given  by: 

N2  +  O2  — ^  2  NO 

as  the  overall  process.  The  elementary  steps  that  determine 
the  rate  of  this  process  depend  on  the  specific  mechanism  of 
NOx  production.  The  relative  importance  in  different  scenar¬ 
ios  of  the  three  popular  mechanisms  of  NOx  production  can 
be  assessed  with  the  implementation  of  such  reduced  chem¬ 
istry  in  the  appropriate  reactive  flow  simulation  models.  As 
an  example,  NO  formation  in  a  Heptane  diffusion  flame  as 
predicted  by  different  models  is  shown  in  Fig.  4.  This  brings 
out  another  valuable  contribution  that  can  be  made  by  com¬ 
putational  combustion.  It  also  shows  that  further  work  needs 
to  be  done  in  improving  the  reduced  chemistry  models. 
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The  systematic  development  of  single-step  or  multi-step 
reduced  chemistry  models  need  a  detailed  reaction  mechanism 
as  the  starting  point.  In  addition,  it  is  important  to  assess  the 
validity  of  the  simplified  models  by  comparing  the  results 
obtained  using  them  to  those  from  detailed  models.  This  is 
often  done  for  a  simplified  unit  problem  such  as  a  laminar 
premixed  or  diffusion  flame  or  a  shock-tube  problem.  There¬ 
fore,  there  is  a  need  to  develop  efficient  ways  to  compute 
multi-species  reactive  flow  problems.  Such  an  approach  is 
discussed  next. 

Parallel  Reactive-flow  Simulations 

As  discussed  previously,  a  reactive  flow  problem  involves 
a  variety  of  physical  and  chemical  processes  that  evolve  at  their 
own  characteristic  time  scales.  An  efficient  approach  to  solve 
such  problems  is  to  use  process  or  time-step  splitting^^  In  this 
approach,  the  various  processes  are  solved  using  appropriate 
algorithms  and  coupled  together.  However,  this  approach  re¬ 
sults  in  the  amount  of  computations  to  be  different  at  different 
grid  points  and  to  vary  with  time.  For  example,  the  chemistry 
solver  will  require  different  number  of  sub-cycles  to  compute 
for  an  overall  time-step  depending  on  the  local  stiffness  of  the 
solution  at  a  specific  grid-point  at  that  time.  Because  of  this 
difference  in  the  amount  of  computations  at  different  grid 
points,  a  simple  domain  decomposition  approach  where  you 
distribute  the  various  grid-points  among  the  processors  is  not 
adequate.  Therefore,  a  new  approach  has  been  developed  in 
collaboration  with  the  University  of  Maryland  [Moon  et  al.^"^]. 

In  this  approach,  the  various  processes  in  the  flame  are 
placed  into  one  of  two  groups: 

Structured,  where  the  computation  is  based  on  structured 
meshes.  Certain  processes,  e.g.  heat  conduction,  are  repre¬ 
sented  by  partial  differential  equations  which  are  discretized 
by  a  finite  volume  technique.  In  these  methods,  the  solution 
procedure  requires  a  good  deal  of  communication  with  neigh¬ 
boring  points,  but  the  amount  of  computation  at  each  point  is 
nearly  the  same. 

Independent,  where  computation  at  each  grid  point  can  be 
carried  out  independently  without  requiring  any  communica¬ 
tion  between  neighboring  points,  e.g.  chemistry.  However,  the 
solution  procedure  may  require  a  different  amount  of  compu¬ 
tation  at  each  grid  point.  This  distribution  of  work  load  may 
change  from  time  step  to  time  step  as  the  flame  evolves. 

In  a  flame  simulation,  the  structured  processes  and  the 
independent  processes  alternate  within  each  time  step.  The 
implication  here  is  that  there  exist  two  conflicting  demands  for 
efficient  parallel  computation  :  (1)  the  necessity  of  regular 
block  partitioning  of  the  grid  to  take  advantage  of  simplicity 
and  efficiency  of  the  structured  computation,  and  (2)  the 
achievement  of  good  load  balance  in  the  independent  proc¬ 
esses.  Two  sets  of  runtime  support  primitives,  Multiblock 


Figure  4. 

Comparison  of  NO  mole  fraction  in  a  Heptane  diffusion  flame 
as  predicted  by  different  chemistry  models.  (Williams.  Ref  12) 
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PARTI  and  CHAOS,  (developed  at  the  University  of  Maryland 
under  an  ONR  initiative  on  Massively  Parallel  Computational 
Science)  have  been  applied  to  cope  with  both  of  the  demands 
successfully.  In  addition,  a  multigrid  method,  based  on  MAD- 
PACK  is  used  as  the  solution  technique  for  the  elliptic  equation 
that  arises  in  the  fluid  convection  model.  To  balance  the  work 
load  of  the  independent  processes,  block-partitioned  data  must 
be  redistributed  across  processors,  and  then  moved  back  into 
the  original  locations  for  the  next  structured  process.  This 
approach  inherently  requires  a  substantial  amount  of  commu¬ 
nication  to  redistribute  data  every  time  step.  New  algorithms 
and  runtime  primitives  which  perform  the  redistribution  of 
work  load  and  minimize  communication  volume  at  low  cost 
have  been  developed  and  tested. 

With  these  advancements,  the  entire  NRL  flame  code 
has  been  parallelized  and  test  problems  computed  using  it 
have  produced  identical  results  to  those  obtained  using  the 
Cray  C-90  version.  For  a  typical  methane  flame  structure 
computation,  it  requires  64  nodes  of  the  Intel  Paragon  or  16 
nodes  of  the  IBM  SP2  to  equal  one  C-90  processor.  Figure 
5  shows  the  scaling  of  the  different  processes  in  terms  of 
the  computational  time  per  time  step  on  different  number  of 
processors.  All  processes  except  the  convection  scale  ex¬ 
tremely  well.  Convection  scales  poorly  because  it  requires 
communication  across  the  entire  grid.  Fluid  convection, 
however,  takes  only  a  small  portion  of  the  overall  computa¬ 
tional  time  and  can  be  minimized  by  considering  larger 
problems  or  problems  involving  more  complex  chemistry 
than  methane.  Thus,  with  the  advances  in  parallel  comput¬ 
ers,  larger  and  more  complex  flame  problems  can  be  solved 
with  such  an  approach. 
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Modeling  Effects  of  Heat  Losses 

Combustion  in  propulsion  systems  occurs  in  confined 
systems  or  at  least  in  systems  with  walls.  Heat  transfer  to  such 
bounding  surfaces  from  the  burning  flame  must  also  be  ac¬ 
counted  for.  Heat  losses  not  only  represent  a  loss  of  useful 
energy  but  may  also  play  a  significant  role  in  determining  the 
detailed  structure,  stability  and  dynamics  of  flames.  Conduc¬ 
tive  heat  losses,  such  as  to  a  flame  holder  or  burner  can  be 
accounted  for  by  the  choice  of  appropriate  boundary  condi¬ 
tions'®.  Information  is  needed  on  the  material  used  for  the 
combustor  and  its  heat  transfer  characteristics.  Surface  chem¬ 
istry  such  as  recombination  of  radicals  at  the  walls  may  also 
need  to  be  included  in  some  studies.  Another  form  of  heat  loss 
is  that  due  to  radiation  from  by-products  of  combustion  such 
as  water,  CO2  and  soot.  It  is  fairly  straight  forward  to  formulate 
a  sub-model  to  account  for  radiative  tosses  if  the  flame  can  be 
assumed  to  be  optically  thin.  This  is  a  useful  approximation 
when  the  self-absorption  of  radiation  in  the  hot  gas  can  be 
neglected.  In  this  formulation  we  need  information  on  the 
Planck  mean  absorption  coefficient  as  a  function  of  tempera¬ 
ture  for  the  species  of  interest.  For  example,  an  eighth-degree 
polynomial  fit  has  been  developed  from  data  available  in  the 
literature  for  the  absorption  coefficient  for  water  and  €62’^. 
Radiation  from  sooting  flames  is  more  complicated  since  the 
optically  thin  assumption  may  not  be  adequate.  More  complex 
radiation  transport  models  have  been  considered  in  the  diffu¬ 
sion  flame  studies  of  Kaplan  et  al.'^.  The  effects  of  using 
different  levels  of  radiation  models  and  their  input  data  are 
currently  being  investigated  under  the  high-energy  fuels  com¬ 
bustion  program**. 

Multi-Phase  Flows 

So  far  the  discussion  has  focused  on  the  combustion  of 
gaseous  fuels.  Although  these  models  are  directly  applicable 
to  pre-vaporized  and  super  critical  systems  in  addition  to 
gaseous  systems,  many  practical  applications  involve  the  com¬ 
bustion  of  liquid  fuel  sprays  and  heterogeneous  gas-particle 
systems.  While  directly  simulating  the  breakup  and  combus¬ 
tion  of  a  multi-component  fuel  jet  injected  into  a  gaseous 
stream  is  beyond  current  capabilities,  significant  progress  has 
been  made  in  the  simulation  of  single  droplet  combustion  and 
the  representation  of  dilute  and  dense  sprays  in  gaseous  back¬ 
ground  flows. 

In  the  computational  studies  of  energetic  fuel  droplet 
combustion,  a  spherically  symmetric  geometry  is  considered 
and  the  solution  domain  is  divided  into  three  regions:  a  liquid 
droplet  core,  a  two-phase  foamy  layer  and  the  surrounding 
gaseous  region”  as  shown  schematically  in  Figure  6.  The 
continuity,  momentum,  energy  and  species  conservation  equa¬ 
tions  are  solved  numerically  in  each  region.  An  exothermic, 
temperature-dependent  reaction  mechanism  is  included  in  the 
liquid  phase  but  mass  diffusion  in  the  liquid  phase  is  neglected. 


The  input  parameters  for  such  a  simulation  include  initial 
droplet  size  and  temperature,  ambient  pressure  and  tempera¬ 
ture,  heat  of  reaction  in  the  liquid,  activation  energy  and 
pre-exponential  constant  for  the  liquid-phase  reaction,  and 
thermophysical  properUes  for  the  gas  and  liquid.  Progress  to 
date  using  such  an  approach  can  be  found  in  Ref.”. 

Information  on  the  vaporization  and  combustion  of  single 
droplets  may  be  used  as  sub-models  in  the  simulation  of  spray 
combustion.  If  the  spray  is  dilute,  it  is  feasible  to  track  the 
motion  of  individual  droplets  or  particles  in  a  flow  system. 
Such  an  approach  has  been  used  to  study  the  dispersion  of 
particles  injected  into  the  high-speed  shear  layer  of  an  axisym- 
metric  combustor.  The  particles’  position  in  the  flow  field  is 
computed  by  solving  the  Lagrangian  equations  of  motion  for 
the  particles  taking  into  account  the  inertial  drag  force  that 
depends  on  the  density  and  size  of  the  particles.  The  flow  field 
is  computed  using  the  MILES  approach  with  the  Flux-cor¬ 
rected  Transport  (FCT)  algorithm.  Details  of  the  model  as  well 
as  more  comprehensive  discussion  of  the  various  cases  simu¬ 
lated  are  documented  elsewhere^**  while  only  highlights  are 
presented  here. 

For  a  typical  case,  an  axisymmetric  jet  with  a  mean 
velocity  of  100  m/s  flows  through  a  cylindrical  inlet  of  diame¬ 
ter  D  into  a  cylindrical  combustion  chamber  of  twice  the 
diameter.  An  annular  exit  at  the  end  of  the  combustion  chamber 
is  modeled  to  produce  choked  flow.  Particles  are  injected  from 
the  inlet-combustor  junction  with  a  stream  wise  velocity  of  50 
m/s  and  zero  radial  velocity.  In  Figure—?,  the  particle  positions 
within  the  combustor  for  particles  of  different  sizes  are  shown 
at  a  particular  time,  20  000  time  steps  (~7 .5  ms)  from  the  start 


Figure  5. 

Scaling  of  different  processes  in  a  methane  flame  calculation 
on  the  Intel  Paragon,  (moon  et  al.,  Ref.  14) 
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of  the  simulations.  This  figure  highlights  the  distinctly  differ¬ 
ent  interactions  between  large-scale  vortical  structures  and 
particles  of  different  sizes.  Very  small  particles  such  as  the  1 
micron  particles  are  engulfed  by  the  vortical  structures  and 
move  with  them.  Very  large  particles  such  as  the  30  micron 
particles  are  not  strongly  affected  by  the  structures  while  the 
intermediate  size  particles  show  a  strong  interaction  in  the 
sense  they  appear  to  form  distinct  sheets  as  they  are  entrained 
along  the  outer  regions  of  the  vortical  structures  and  then  flung 
out  These  observation  support  the  different  mechanisms  for 
particle  dispersion  postulated  in  the  earlier  study  of  particle 
dispersion  in  axisymmetric  jets^^  Another  significance  of  this 
study  is  on  the  choice  of  particle  sizes  used  in  flow  diagnostics 
to  track  the  fluid  flow.  It  is  also  relevant  to  the  combustion  of 
multi-component  fuels  with  additives  that  may  remain  in  solid 
form  in  a  gaseous  background.  For  the  combustion  of  dense 
sprays,  the  interactions  among  the  different  particles/droplets 
must  also  be  taken  into  account.  Furthermore,  it  is  not  compu¬ 
tationally  practical  to  follow  individual  particles  when  there 
are  a  very  large  number  of  them  and  they  are  of  different  sizes. 
The  traditional  method  of  modeling  polydisperse  sprays  is  to 
partition  the  initial  droplet  size  distribution  into  bins  and  fix 
these  bins  so  that  during  the  calculation  as  the  droplets  evapo¬ 
rate,  they  can  change  bins.  This  becomes  difficult  for  dense 
sprays  because  at  a  given  location,  the  range  in  droplet  sizes 
can  be  very  large.  A  variation  of  this  approach  is  currently 
being  considered^^.  In  this  approach,  only  the  initial  droplet 
size  distribution  is  partitioned  into  finite-size  bins.  Each  size 
is  then  followed  through  the  conservation  equations  in  its  own 
system  of  coordinates  which  moves  with  the  droplets.  Thus 


Figure  6. 

Schematic  of  a  vaporizing  liquid-fuel  energetic  droplet  shows  (a) 
the  initial  droplet  at  fi  =  0,  while  (b)  and  (c)  show  the  droplet  at 
progressively  longer  times,  tzandb-  (Sirignano  et  al,  Ref  19) 
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Figure  7. 

Instantaneous  spatial  distribution  of  vorticity  and  particles  of 
different  sizes  in  an  axisymmetric  ramjet  combustor. 
(Kailasanath  etal.,  Ref  18) 
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each  initially  mono-size  bin  develops  a  continuum  of  droplet 
sizes  as  the  droplets  evaporate.  Parametric  studies  using  this 
model  suggest  that  polydispersity  is  beneficial  to  spray  com¬ 
bustion  in  that  the  smaller  droplets  provide  the  fuel  vapor  to 
initiate  the  flame  whereas  the  large  droplets  provide  the  means 
for  the  fuel  to  penetrate  into  the  oxidizer  and  thus  propagate 
the  flame^^. 

The  progress  in  the  modeling  of  particles  and  droplets  in 
fluid  flows  will  also  be  useful  in  the  simulation  of  the  forma¬ 
tion  and  destruction  of  soot. 

Combustion  of  Sooting  Fuels 

Most  often  combustion  in  propulsion  systems  takes  place  in 
a  diffusion  flame  scenario  and  most  of  the  hydrocarbon  fuels  in 
use  in  the  Navy  have  a  tendency  to  soot.  Particularly  so  are  the 
new  class  of  high-energy  fuels  being  developed  which  have  a 
large  carbon  to  hydrogen  ratio,  and  hence  an  increased  sooting 
tendency.  Reduction  of  soot  is  crucial  to  reduce  signatures  from 
propulsion  system  exhausts  as  well  as  to  improve  their  efficiency. 
Numerical  simulation  of  sooting  flames  is  still  in  its  infancy. 
Better  chemical-mechanistic  understanding  of  the  formation  and 
destruction  of  soot^  is  needed  for  the  development  of  sub-models 
for  inclusion  in  numerical  simulations. 
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As  a  first  step,  sooting  methane-air  flames  have  been 
simulated  using  a  flame  sheet  approximation^.  For  these 
simulations,  additional  equations  for  the  transport  of  soot 
volume  fraction  and  number  density  have  been  formulated. 
The  soot  volume  fraction  is  computed  as  a  function  of  the  local 
gas  properties  (temperature,  density,  fuel  mole  fraction)  using 
a  simplified  model  that  is  based  on  experimental  measure¬ 
ments^^.  However,  the  model  does  consider  the  effects  of 
nucleation,  surface  growth,  coagulation,  and  oxidation. 

In  spite  of  the  simplifications  adopted  in  the  model,  com¬ 
parisons  between  the  computations  and  experiments  show 
reasonable  agreement  in  the  overall  location  of  the  flame 
surface  and  the  flame  height.  The  temperature  and  soot  volume 
fractions  are  also  in  reasonable  agreement  at  most  measure¬ 
ment  locations.  Although  the  comparisons  are  encouraging, 
further  studies  need  to  be  made  in  more  complex  combustion 
scenarios. 

Concluding  Remarks 

In  this  article,  we  have  highlighted  the  unique  computa¬ 
tional  issues  that  need  to  be  addressed  to  simulate  combustion 
problems  relevant  to  the  Navy’s  propulsion  program.  We  have 
briefly  reviewed  the  recent  progress  focusing  on  the  work 
sponsored  by  the  Navy,  either  directly  by  ONR  or  indirectly 
through  NRL.  The  work  takes  advantage  of  developments  in 
CFD  but  as  discussed  above  there  are  also  several  unique 
aspects  of  combustion  that  must  be  addressed  directly.  Ad¬ 
vances  in  computers  and  progress  in  computational  methodol¬ 
ogy  will  continue  to  expand  the  range  of  problems  that  can  be 
tackled.  Because  of  the  complexity  of  the  combustion  phe¬ 
nomenon,  modeling  of  turbulence  and  other  processes  such  as 
soot  formation  and  destruction,  radiation  and  chemistry  will 
be  required.  With  the  increased  emphasis  on  timed  fuel  injec¬ 
tion  for  increased  fuel  efficiency  and  reduced  pollution,  more 
emphasis  will  be  placed  on  modeling  time-dependent  turbu¬ 
lent  combustion.  Meanwhile,  simulation  of  more  idealized 
systems  and  scenarios  have  already  made  an  impact  in  improv¬ 
ing  our  current  understanding  of  the  dynamics  and  control  of 
combustion.  This  is  also  paving  the  way  for  optimized  design 
of  complex  combustion  systems. 
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